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ABSTRACT 

The  equations  for  estimating  the  performance  of  an  axial  flow 
turbine  are  examined  to  determine  the  efficiency  and  other  parameters 
when  very  high  power  outputs  per  stage  are  required.   It  is  indicated 
that  turbine  stages  should  be  designed  for  peak  efficiency  or  for 
high  power  output,  but  that  compromise  designs  are  not  advisable. 

A  survey  of  the  immediately  available  literature  on  losses  in 
turbine  stages  at  the  U.  S.  Naval  Postgraduate  School  is  outlined. 
An  attempt  to  find  and  correlate  specific  loss  data  with  theoretical 
relations  was  not  successful.   Data  is  fragmentary,  theoretical  re- 
lations are  few,  and  the  correlation  of  the  available  data  on  profile 
losses,  tip  clearance  losses,  effect  of  trailing  edge  thickness  and 
secondary  losses  is  poor. 

An  axial  flow  single  stage  turbine  was  designed  for  transonic 
flows  and  a  digital  computer  program  was  written  to  evaluate  the 
off-design  performance  of  this  turbine.   The  particular  relations  for 
supersonic  expansion  after  a  blade  row  are  included  in  the  program 
as  well  as  the  effects  of  entry  shocks  into  the  rotor.   Initial 
indications  are  that  the  program  should  be  useful  for  evaluating 
turbines  designed  by  the  procedure  shown. 

The  computer  programs  included  are  written  in  FORTRAN  language 
for  the  CDC  1604  Digital  Computer. 
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THE  PERFORMANCE  OF  HIGH  POWER  OUTPUT  AXIAL 
FLOW  TURBINES  UTILIZING  FLOWS  IN  THE  TRANSONIC  REGIME 

I.   Introduction 

The  increasing  need  in  space  applications  for  turbines  of  high 
power  output  from  a  minimum  physical  volume  and  weight  requires  more 
knowledge  of  efficiency,  losses  and  performance  estimation  of  these 
turbines. 

If  the  required  power  output  and  space  limitations  are  predominant 
then  efficiency  must  suffer.   As  a  means  of  showing  the  degree  of 
efficiency  sacrificed,  the  equations  which  determine  the  power  output 
and  efficiency  of  an  arbitrary  turbine  were  programmed  for  digital 
computer  solution  in  terms  of  the  head  coefficient.   The  head  coef- 
ficient for  peak  efficiency  was  determined  and  then  the  head  coef- 
ficient was  increased  in  increments  and  various  parameters  affecting 
the  turbine  performance  were  computed  at  each  increment  until  the 
efficiency  dropped  to  an  arbitrary  cutoff  point.   The  data  are  plotted 
in  non-dimensional  form  and  simple  equations  are  found  to  represent 
the  data. 

Of  paramount  importance  in  the  design  of  a  turbine  is  the  esti- 
n^ition  of  the  losses  due  to  boundary  layer  formation,  mixing  effects, 
clearances  and  secondary  flows.   This  area  is  explored  in  the  hope 
of  finding  a  more  accurate  and  precise  method  of  presenting  these 
losses.   A  survey  of  the  immediately  available  literature  is  included. 


A  particular  single-stage  turbine  is  designed  for  supersonic 
flow  leaving  the  stator  blade  row  and  a  relative  subsonic  flow  as 
seen  by  the  rotor.   The  equations  for  determining  the  effects  of 
after-expansion  behind  a  blade  row  and  entry  shocks  entering  a 
blade  row  are  programmed  for  the  computer  in  addition  to  the 
regular  performance  parameters.   The  rpm  and  pressure  ratio  across 
the  turbine  are  then  varied  to  show  the  effects  of  demanding  extreme 
power  outputs  from  a  turbine  so  designed.   A  performance  map  of  the 
turbine  is  developed. 
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II.  Development  of  High  Head  Coefficient  Data 

A.   Basic  Equations 

In  Ref.  1  the  equations  are  developed  to  establish  the  design 
parameters  for  turbine  stages.   These  equations  are  obtained  by 
assuming  that  the  flow  through  the  turbine  is  represented  by  the 
conditions  on  a  mean  stream  surface  and  the  primary  parameter  is 
taken  to  be  the  "head  coefficient".   The  head  coefficient  relates 
the  total  isentropic  drop  of  the  stage  from  the  inlet  stagnation 
pressure  to  the  discharge  static  pressure  to  the  kinetic  "energy" 
of  the  peripheral  speed  at  the  rotor  entrance  and  is  defined  as 

J  _    £^=    =  f^\ 
-*"-    tj*/2        [a  <» 


U*/2 


In  order  to  relate  the  theoretical  expressions  to  values  found 
in  engineering  use,  some  expression  for  the  losses  involved  must  be 
used.   A  velocity  coefficient  (  )r   )  is  used  as  a  measure  of  rotor 
efficiency  and  is  defined  as 

V/z 


WziL  (2) 


This  velocity  coefficient  has  been  found  through  experience  to  depend 
primarily  on  the  deflection  of  the  flow  in  the  rotor  {£jo   )■   Vavra 
has  established  a  mean  curve  through  data  from  a  variety  of  sources 


r  /o4  i        /80-A/8 


/80-AjS 

(3) 
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and  this  equation  is  used  to  represent  the  losses  for  purposes  of 
evaluating  the  efficiency  at  each  step  of  the  calculation. 

Following  the  procedure  outlined  in  Ref.  1,  arbitrary  values 

were  chosen  for  the  angle  of  the  flow  leaving  the  stator  (Of,), 

reheat  factor  (f),  and  the  carryover  coefficient  of  kinetic  energy 

into  the  next  stage  (©  ),  the  radius  ratio  (R2/R-) ,  the  meridional 

velocity  ratio  (V  0/V  ,),  the  velocity  coefficient  in  the  stator  (r), 
m/  ml  ' 

and  the  carryover  coefficient  of  kinetic  energy  from  the  stator  into 

the  rotor  (&„)•      Then  the  degree  of  reaction  (r  )  is  selected  and 

the  head  coefficient  is  increased  in  increments.   The  degree  of 

reaction  is  the  fraction  of  the  isentropic  enthalpy  change  through 

the  turbine  stage  which  occurs  in  the  rotor.   At  each  increment  of 

head  coefficient  the  particular  values  of  velocity  coefficient  (j), 

internal  efficiency  (/j)i  leaving  loss  coefficient  (k_),  ratio  of 

meridional  absolute  velocity  to  peripheral  speed  (V  , /Un ) ,  relative 

ml   1 

velocity  ratio  (W^/W.),  and  the  relative  and  absolute  angles  (0^2> 


k? 


_)  are  computed. 


A  similar  computer  program  using  the  same  equations  was  written 
to  find  the  peak  efficiency  for  a  given  entering  angle  and  degree  of 
reaction,  and  then  all  the  other  parameters  are  computed.   The  values 
are  also  found  for  higher  head  coefficients  such  that  the  efficiency 
is  0.57„  less  than  optimum  and  then  1.0%  less  than  optimum. 
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Computer  Program 
The  equations  to  be  solved  by  iteration  until  y^    my   are; 

/ 


1.         AAyy^Si    =    ^Osy*sO(i  


Mlw  Pz   


•sUls 


f  (%*/&>;)    ax-0*/  "/TvJ*1 


<*) 


(5) 


2       -i 


where 


Xz=  ^^v^[4^^6^Ja*» 


-2<fo  ^Mn-Of,  J 


/-A*  I 


**$ 


/ 


(5a) 


3-  #     =    0#- 


^.^<9 


y- 


4.97 


'a  /o4   r      /80   Aft 

Once  this  condition  (  j^    ■  jr  )  is  satisfied  the  values  of  the 
other  parameters  are  found  from 


(3) 


*w2 


(6) 


(7) 


& 


"^ 


/    /A1* 


/X*      U/      ^K 


=  r^-JE  <<***, 


(8) 
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The  equations  were  coded  in  Fortran  language,  Ref.  2,  and  solved 
on  the  Control  Data  Corporation  1604  computer  at  the  U.  S.  Naval 
Postgraduate  School.   A  copy  of  the  program  is  presented  in  Appendix  III. 
The  resulting  data  printouts  are  included  in  Appendix  V  as  Table  I. 

The  program  finding  the  peak  efficiencies,  efficiencies  0.57«  less 
than  optimum  and  efficiencies  1.07„  less  than  optimum  is  also  in 
Appendix  III,  with  data  in  Appendix  V  as  Table  II.   It  should  be 
noted  here  that  all  equations  are  completely  coded  in  Fortran  language 
so  that  any  selected  numerical  values  may  be  introduced  in  the  be- 
ginning of  the  program.   Any  of  these  arbitrarily  chosen  values  may 
be  changed  and  the  equations  of  the  main  program  remain  unchanged. 

C.   Results 

Representative  data  from  Tables  I  and  II  are  presented  in  Fig.  1. 
The  regularity  and  progression  of  the  data  curves  indicate  that  non- 
dimensionalizing  the  efficiency  and  head  coefficient  parameters  might 
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I e  profitable.   When  all  the  efficiencies  for  one  entering  angle  and 
one  degree  of  reaction  are  divided  by  the  optimum  efficiency  for  that 
<^et  of  conditions  (  yf   /ODt)  and  the  head  coefficients  are  divided 

by  the  optimum  head  coefficient  (k.  /k.     )  the  curves  of  Figs.  2, 

xs   xs  op L 

3,  4  and  5  result. 

A  closer  inspection  of  these  plots  reveals  that  for  each  degree 
of  reaction,  one  curve  suffices  for  all  four  entering  angles. 
The  cubic  equation 

tf'JopJ   -  -000223(^ ]     +.008345(^ J     -.uow/JiS )+1.125 

I    is  °P'/  1     is  opt/  1  kls  opt 


it 

represents   all   the   data   for   r     *  0.0   for  all   entering  angles    from  60° 

to   75°. 


Similarly  as   an  approximation   to  all   the   data  at   r     =0.25,    the 


equation   is 


(7/7opt>   '  -000223  [J±2 J'     +.00770      [{H ]'     -  .  1098   /  £* 


-r.uu/zu -.1098       ~= +l.r 

is  opt/  \     is   opt  /  Kis  opt 


and   for   r     *  0. 50 


'7/7oPt> 


-.000125 


kis  3  /k^  \  2  /k- 

+.00554  MS -.09414  MS 

\    is  °Pt/  I     is  opt  I  I  Kis  opt 


The  maximum  error  in  the  range  of  entering  angles  presented  here  is 
approximately  8  parts  in  130.  This  is  probably  accurate  enough  for 
preliminary  engineering  estimates   based  on   the   equations    in  Ref.    1 
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used  in  the  computer  programs.   As  will  be  seen  later,  the  loss  data 
is  not  known  to  this  accuracy. 

The  curves  show  that,  roughly  speaking,  increasing  the  head 
coefficient  to  five  times  the  optimum  causes  a  reduction  of  307<>  of 
the  obtainable  peak  efficiency.   Increasing  the  head  coefficient  to 
ten  times  the  optimum  causes  a  reduction  of  457o  of  the  obtainable  peak 
efficiency.   This  indicates  that  if  it  is  necessary  to  go  to  very  high 
head  coefficients  to  obtain  a  certain  required  work  output,  the  rate 
of  change  of  efficiency  with  increasing  head  coefficient  is  greatest 
from  optimum  to  five  times  optimum  head  coefficient.   The  rate  of 
change  of  efficiency  decreases  as  the  head  coefficient  is  further 
increased. 

Table  I  of  Appendix  V  can  be  used  to  establish  the  other  stage 
parameters  after  the  head  coefficient  range  of  interest  is  chosen. 

Figs.  1  through  5  are  not  the  off-design  performance  of  a 
particular  turbine,  but  they  are  the  locus  of  points  of  an  infinite 
number  of  different  turbine  designs.   Each  point  on  each  curve 
represents  a  turbine  designed  for  that  condition.   Therefore,  if  the 
requirement  is  that  the  turbine  be  designed  for  a  high  head  coef- 
ficient, the  loss  in  efficiency  need  not  be  the  absolute  controlling 
factor,  and  the  design  may  as  well  be  based  on  extreme  high  values 
of  head  coefficient.   These  investigations  show  that  it  is  advisable 
to  either  design  the  turbine  for  peak  efficiency  of  maximum  output 
but  do  not  compromise  both  by  operating  at  only  two  or  three  times 
the  optimum  head  coefficient. 
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III.   Loss  Investigation 

Since  the  design  of  a  turbine  depends  on  the  loss  assumptions 
made,  the  use  of  an  average  velocity  coefficient  based  on  general 
considerations  of  flow  deflection  in  the  rotor  is  not  a  very  precise 
method  of  estimating  the  losses.   This  was  the  method  used  in  the 
preceding  section  of  this  thesis.   A  better  way  should  be  found  and 
the  following  considerations  are  basic  to  the  problem. 

A.   Loss  Components 

It  seems  logical  to  class  the  losses  encountered  in  turbines 
into  four  particular  categories; 

1)  Profile  loss  --  the  friction  loss  associated  with  the 
rcrmation  of  the  boundary  layer  on  the  blades  of  the  cascade. 

2)  Tip  clearance  loss  --  the  pressure  and  energy  loss  associated 
with  the  leakage  flow  between  the  tips  of  the  blades  and  the  casing 

of  the  rotor,  including  the  effects  of  trailing  vortices  at  the  hub 
rtiere  there  is  zero  clearance. 

3)  Trailing  edge  loss  --  the  loss  due  to  mixing  effects  behind 
the  blades  of  the  row  as  caused  by  the  finite  thickness  of  the  blade 
trailing  edge. 

4)  Secondary  loss  --  the  loss  due  to  the  vortices  which  are  set 
up  in  the  flow  in  the  curved  channel  between  blades  giving  a  velocity 

^onent  normal  to  the  desired  direction. 
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In  the  axial  flow  turbine,  the  closely  packed  blades  and  large 
flow  deflections  usually  require  the  flow  through  the  turbine  to  be 
considered  from  the  point  of  view  of  "channel  flow"  between  blades 
instead  of  using  airfoil  theories.   As  will  be  seen  in  the  subsequent 
discussion,  this  secondary  flow  causes  the  principal  problem  pre- 
venting the  accurate  assessment  of  the  losses. 


Looking  at  the  flow  between  two  blades  of  an  arbitrary  cascade, 
Fig.  6,  the  velocity  distribution  is  as  shown  for  the  mid-span  of 
the  blade  channel  and  it  can  be  seen  that  a  static  pressure  gradient 
is  developed  between  surfaces  0-1-2  and  0,-l'-2l.   The  centrifugal 
forces  on  the  fluid  being  turned  cause  an  increase  in  the  static 
pressure  near  the  0'-l'-2'  wall  and  decrease  the  static  pressure 
near  the  0-1-2  wall. 


/#•  J7 

'      *' 

/     J7 

/      */ 

Fig.  6  Flow  Between  Blades  of  an  Arbitrary  Cascade 
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Assuming  non-viscous  flow  outside  the  boundary  layer  requires  the 

flow  to  increase  velocity  near  the  0-1-2  convex  wall  to  maintain  the 

9 
total  pressure  a  constant  (P  "  p  +  %  Ow" ).   Conversely,  the  flow 

velocity  must  decrease  near  the  O'-l'-Z'  concave  wall  outside  the 

boundary  layer.   In  the  boundary  layer  at  the  top  and  bottom  of  the 

channel,  the  loss  of  velocity  (and  total  pressure)  due  to  friction 

causes  the  flow  direction  to  be  dominated  by  the  pressure  gradient 

imposed  (assuming  the  static  pressure  is  constant  through  the  boundary 

layer).   Therefore,  as  the  flow  progresses  from  station  0-0'  to 

station  1-1'  to  station  2-2'  ,  the  fluid  in  the  upper  and  lower 

boundary  layers  tends  to  flow  from  0'  toward  1  and  from  1'  toward  2. 

This  change  of  direction  of  the  flow  in  the  boundary  layers 
results  in  the  formation  of  two  vortices  superimposed  on  the  main 
flow  leaving  the  channel  at  station  2-2'.   In  Fig.  7  the  flow  patterns 
in  the  upper  and  lower  boundary  layers  are  shown  as  well  as  the 
induced  vortices  rotation  direction  at  the  exit  plane. 

This  vorticity  results  in  velocity  components  developed  perpen- 
dicular to  the  desired  exit  velocity.   Since  a  velocity  component 
normal  to  the  desired  direction  cannot  be  recovered  usually,  most  of 
the  useful  energy  in  these  components  is  lost  in  friction.   This  is 
the  so-called  "secondary  loss".   For  a  rotating  row  of  blades,  the 
blade  tip  clearance  also  has  a  marked  effect  on  the  secondary  flow. 
The  leakage  from  the  high  pressure  to  the  low  pressure  side  of  a  blade 
encourages  an  additional  trailing  vortex  to  form  a  reaction  blading. 
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For  impulse  blading  the  flow  over  the  tips  which  is  not  deflected  also 
enhances  vortex  formation.   For  cascades  of  infinite  blade  height 
(two-dimensional  only)  there  is  no  secondary  flow  since  the  mechanism 
depends  principally  on  the  annulus  boundary  layers. 


Fig.  7  Flow  in  the  Boundary  Layer  of  a  Cascade  Channel 

With  the  physical  parameters  of  blade  height,  trailing  edge 
thickness,  tip  clearance,  profile  shape,  and  spacing  given,  it  should 
be  possible  to  make  a  reasonably  exact  estimate  of  losses. 


B.   Survey  of  the  Literature 

The  Immediately  available  literature  at  the  U.  S.  Naval  Post- 
graduate School  was  surveyed  to  find  the  particular  loss  breakdown 


used  by  the  various  sources,  the  theory  and  representations  offered, 
and  any  precise  loss  data  published.   The  survey  is  limited  in  scope 
and  critical  since  the  particular  desire  is  to  locate  specific  loss 
theories  and  data.   The  results  of  this  survey  are  presented  here 
in  the  order  of  reference  books,  general  texts,  and  then  periodicals 
and  reports.   The  NACA/NASA  reports  provide  a  wealth  of  data  on  the 
efficiency,  or  total  loss  coefficient,  of  specific  turbine  designs 
from  tests.   Very  little  information  was  found  for  loss  components, 
however,  in  the  frame  of  reference  considered  here. 

1.   In  Ref.  1  an  equation  is  derived  for  the  determination  of 
a  flow  coefficient  in  a  rotor  i^r)   as  a  function  of  the  change  in 
flow  angle  of  the  relative  flow.   This  equation  is  a  fitted  curve 
through  data  from  several  sources  and  a  variety  of  blade  shapes  and 
operating  conditions  and  is  the  same  equation  used  in  the  preceding 
section  of  this  thesis.   In  Fig.  8  is  shown  the  general  shape  of  the 
curve  and  the  derived  equation  (3). 

A  loss  coefficient  can  be  defined  as 

,Z 


f-/-r 


(id 


Considering  the  case  of  A  a    ■  0,  we  see  the  loss  coefficient 
7   ■  1  -  j  represents  those  losses  in  the  cascade  which  are 

present  when  the  flow  is  not  deflected.   As  the  flow  is  deflected 
an  amount  &p   ,  the  reduction  in  ^indicates  an  additional  loss. 
The  additional  loss  will  be  taken  as  the  secondary  loss.   The 
secondary  loss  coefficient  is  then 

*         ,Ul  (12) 


i=f-?0  =i-f 
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Curves  shown  in  subsequent  discussions  of  secondary  losses  will  now 
show  a  loss  coefficient  curve  based  on  this  equation  plotted  as  a 
dashed  line  and  labeled  "VAVRA".   This  curve  will  act  as  a  general 
reference  and  will  show  Maverage"  conditions  so  a  comparison  of  each 
method  can  be  made. 


/.a 
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Fig.  8  Flow  Coefficient  r  as  a  Function 
of  Turning  Angle  in  a  Cascade 


2.   A  veritable  fountainhead  of  theory,  design  and  ideas  are 
presented  in  Ref.  3.   This  basic  work  by  Stodola  touches  on  almost 
all  the  subjects  considered  important  now.   The  important  works  of 
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others  up  until  1927  are  included  to  show  their  measurements  and  de- 
velopments.  The  losses  through  fixed  rows  of  blades  are  given  in 
the  form  of  the  so-called  velocity  coefficient  j    \    the  ratio  of  actual 
velocity  to  theoretical  velocity.   The  measurements  of  the  Institution 
of  Mechanical  Engineers  of  England  and  others  all  show  good  correlation 
for  flows  up  to  sonic  speed  for  stators  with  an  axial  entering  angle 
and  an  exiting  angle  of  70  -78   from  the  axial.   For  supersonic 
speeds  the  data  show  some  scatter  due  to  difference  in  measuring 
apparatus  and  some  basic  differences  in  design  of  the  convergent- 
divergent  passage. 

Secondary  flows  in  bends  are  explained  but  there  is  wide  vari- 
ation in  the  available  data  for  friction  losses  in  bends.   No  pub- 
lished data  were  available  for  the  losses  in  rotating  blade  rows. 
The  data  and  experiments  on  fixed  rows  of  impulse  rotor  blades  were 
of  limited  application  and  only  show  general  trends. 

For  design  purposes,  velocity  coefficients  and  loss  coefficients 

y 

are  given  in  general  terms  only,  i.e.,  "velocity  coefficient  for 
stators  ■  0.975-0.92  for  long  and  short  nozzles  respectively".   The 
curve  presented  by  Stodola  for  the  velocity  coefficient  for  rotating  , 
blade  rows  as  a  function  of  turning  angle  is  included  in  the  data 
used  for  the  determination  of  the  equation  (3)  by  Vavra . 

The  tip  clearance  losses  had  been  explored  in  a  limited  fashion. 
These  losses  were  accounted  for  by  an  empirical  formula  for  the  de- 
crease in  internal  work  as  a  function  of  clearance  and  blade  height, 


15- 


• 


6.27  k   /h  (13) 

where  k  "  radial  clearance  and  h  is  the  blade  height  corresponding 
to  zero  clearance.   This  relation,  credited  to  Anderhub,  is 
plotted  on  Fig.  9  and  noted  as  (A) . 

The  steam  turbine  design  procedures  presented  are  strongly 

dependent  on  making  the  new  design  a  small  modification  of  an  already 

designed  unit  which  is  in  service  and  for  which  the  efficiency 

is  known.   The  peak  efficiencies  were  recognized  to  be  in  the  regions 

of  k,   from  1.0  to  5.0,  as  shown  in  the  previous  section  of  this 
is  r 

paper. 

3.  The  axial  turbine  section  of  Ref.  4  gives  a  lucid  and  brief 
explanation  of  the  flows  in  such  machines  but  the  loss  data  and 
secondary  flow  considerations  are  based  completely  on  Ref.  22,  23 
and  24  which  are  considered  later  in  this  paper. 

4.  Since  books  written  to  be  used  as  school  texts  of  necessity 
try  to  cover  the  complete  spectrum  of  fluid  mechanics,  theory,  and 
design,  only  the  most  cursory  reference  is  made  to  specific  losses 
in  turbine  stages.   Refs.  5,  6  and  7  are  representative  of  books  of 
this  type.   Ref.  5  defines  a  work  recovery  coefficient  as  the  ratio 
of  useful  work  extracted  by  the  rotor  to  the  total  kinetic  energy 

of  the  flow  leaving  the  stator,  times  the  stator  flow  efficiency. 
This  work-recovery  coefficient  variation  with  the  velocity  ratio 
(U/V..)  is  shown  for  impulse  and  reaction  conditions  and  no  other 
explanation  of  the  nature  of  the  losses  is  offered.   The  charts 
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resented  are  based  on  data  from  Ref.  16,  in  which  report  the  losses 
in  a  blade  row  are  presented  only  as  a  function  of  stagger  angle  and 

incidence. 

In  Ref.  6  the  author  makes  the  statement  that  no  adequate 
criterion  is  available  for  the  difference  in  the  flow  conditions 
between  an  efficient  and  an  inefficient  hydrodynamic  machine  of  a 
given  type.   Any  losses  and  efficiencies  of  turbines  offered  are 
tied  to  airfoil  theory  or  definitions  of  losses  by  Stodola  and 
others.   Secondary  flows  and  radial  flow  of  the  boundary  layer  are 
mentioned  only  from  Ref.  17  where  the  measurements  were  made  on 
widely  spaced  blades  on  a  rotating  blade  row  at  low  speed. 

A  large  list  of  references  on  the  work  in  gas  turbines  is 
given  in  Ref.  7   for  the  data  available  at  that  time.   Principally 
the  work  in  Ref.  16  below  is  mentioned  in  connection  with  turbines 
but  only  general  terms  are  offered  and  no  specific  loss  derivations 
are  offered. 

5.   The  most  comprehensive  effort  to  make  a  breakdown  of  losses 
and  show  specific  data  found  in  a  textbook  is  that  of  Ref.  8.   This 
book  takes  selected  data  from  Refs.  15,  16  and  23,  discussed  later 
in  this  paper,  and  presents  plots  of  the  losses  due  to  leakage, 
profile  losses,  losses  from  turbulence  and  wall  friction  (including 
secondary  loss  effects)  and  incidence  losses.   The  discussion  is 
brief,  however,  and  no  new  data  or  theory  are  presented. 
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6.  An  attempt  to  mathematically  tackle  the  problem  of  the 
secondary  flows  in  bends  is  presented  in  Ref.  9.   The  equations  of 
motion  of  an  ideal  fluid  are  used  and  frictional  effects  are  ignored. 
An  induced  drag  pressure  drop  coefficient  is  defined  to  account  for 
the  power  loss  due  to  secondary  flow.   The  method  is  applied  to  the 
flow  in  the  corner  of  a  3'  x  4'  wind  tunnel  which  has  two  turning 
vanes  to  simulate  a  blade  row  and  the  method  shows  the  proper  trends. 
A  quantitative  value  derived  for  the  loss  is  computed  and  found  to 
be  about  twice  the  measured  loss.   No  other  example  of  the  theory 
being  applied  is  available. 

7.  A  qualitative  investigation  of  the  flow  in  boundary  layers 
and  wakes  of  blade  rows  is  offered  in  Ref.  10.   An  excellent  discussion 
of  the  general  nature  of  these  flows  is  accompanied  by  the  statement 
that  a  quantitative  theory  for  these  effects  is  not  available. 

8.  An  analytical  method  of  estimating  turbine  performance 
described  in  Ref.  11  is  based  on  losses  due  to  incidence  angle  and 
losses  due  to  flow  of  the  fluid  through  the  blade  row.   A  blading 
loss  parameter  is  used  to  represent  the  losses  but  only  a  general 
range  of  values  of  the  parameter  are  given.   No  method  is  presented 
for  determining  the  magnitude  of  the  blading  loss  parameter  before 
the  turbine  is  built. 

9.  The  general  state  of  the  art  to  1948  is  described  in  Ref.  12. 
The  paper  is  not  intended  to  be  used  as  a  basis  for  other  than  general 
design  considerations,  however,  and  some  empirical  relations  presented 
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were  attempts  to  cover  both  axial  turbines  and  compressors.   This 
necessarily  limited  the  worth  of  such  relations  for  large  flow 
deflection  turbine  loss  estimations.   Those  relations  presented  for 
strictly  turbine  use  are  repeated  in  Refs.  24  and  25  which  will  be 
covered  in  more  detail  later  in  this  thesis. 

10.  The  losses  in  stators  are,  of  course,  easier  to  measure 
than  the  losses  in  rotors.  Some  excellent  tests  and  measurements 
are  presented  in  Refs.  13  and  14.  These  show  that,  in  general,  a 
properly  designed  nozzle  for  a  turbine  will  have  a  velocity  coef- 
ficient (ratio  of  actual  velocity  to  isentropic  velocity)  value  of  .96, 
This  single  value  is  found  to  be  uniform  over  a  considerable  range 

of  designs  and  conditions.   Off-design  values  do  fluctuate,  partic- 
ularly in  the  case  of  supersonic  nozzles  but  the  reasons  do  not  seem 
obscure  and  loss  estimations  for  stator  rows  should  be  on  relatively 
firm  ground. 

11.  Ref.  15  is  mentioned  as  a  reference  work  in  many  textbooks 
on  turbomachinery.   This  compilation  of  the  then  current  design 
practices  only  shows  the  losses  as  a  function  of  axial  length  of  the 
blades  of  the  rotor,  turning  angle  in  the  rotor  and  incidence  angle. 
The  curves  for  loss  as  a  function  of  turning  angle  are  included  in 
the  data  used  for  the  determination  of  the  equation  (3)  by  Vavra. 
Specific  detailed  loss  breakdowns  or  data  are  not  available. 

12.  The  systematic  investigations  into  secondary  loss  presented 
in  Refs.  18  and  19  are  of  real  interest  although  the  test  data  is  too 


19- 


limited  to  permit  any  universal  application.   An  attempt  is  made  to 
define  the  two-dimensional  loss  by  theory,  measure  the  end  loss  at 
the  blade  root  and  the  clearance  loss  at  the  tip,  and  subtract  these 
from  the  total  pressure  loss  measured  to  arrive  at  the  secondary 
loss.   The  tests  are  performed  on  a  particular  airfoil  shape  in 
three  solidity  arrangements  and  three  flow  deflections.   The  inac- 
curacies and  errors  discussed  should  help  any  further  investigator 
in  this  subject. 

13.   A  comparable  breakdown  of  the  losses  into  components  as 
assumed  here  is  offered  in  Ref.  20.   This  work  by  Markov  also  includes 
the  experimental  results  of  some  other  Russians  in  the  field.   With- 
out attempting  to  completely  describe  the  development  of  equations 
and  data  presented,  the  pertinent  data  for  the  present  consideration 
is  extracted. 

The  profile  loss,  or  two-dimensional  loss  which  would  be  caused 
by  blades  of  infinite  length  is  shown  on  Fig.  10  along  with  other 
data  to  be  described  later.   The  basis  of  the  Ref.  20  data  is  a  whole 
series  of  experimental  measurements  of  widely  used  cascades  for 
turbines  in  Russia.   For  stators,  the  profile  loss  coefficient  is 
almost  constant  at  .03  over  a  wide  range  of  conditions. 

The  tip  clearance  loss  is  shown  on  Fig.  9  as  represented  by  two 
plots  in  Ref.  20.   One  plot  is  shown  as  the  loss  for  a  reaction  stage 
and  one  plot  is  based  on  experimental  data  for  the  efficiency  of  stages 
with  a  small  degree  of  reaction.   The  magnitudes  of  the  degree  of 
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reaction  are  not  specified  but  the  trend  is  correct.   Stages  with  a 
high  degree  of  reaction  have  a  higher  loss  across  the  tip,  as  would 
be  expected  due  to  the  pressure  difference  across  the  rotor  from 
entrance  to  exit. 

The  trailing  edge  loss  for  cases  where  the  pressure  gradient 
is  small  along  the  surface  of  the  trailing  edges  of  the  blades  is 
computed  by 


(14) 


Where  such  a  pressure  gradient  is  present  the  use  of  boundary  layer 
thickness  parameters  is  required.   The  trailing  edge  loss  coefficient 
according  to  eq.  (14)  is  shown  in  Fig.  20. 

The  energy  lost  due  to  the  secondary  flow  is  treated  as  a 
function  of  the  mass  flow  rate  in  the  cascade.   The  efficiency  of  a 
cascade  of  blades  with  finite  length  is  defined  as 

*V  m  4y  .  Energy  lost  due  to  secondary  flow I 

/(2-dim)        Energy  available  flowing  through  cascade! 

(15) 
The  loss  coefficient  is  defined  as 


r=/-; 
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and  so 
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The  secondary  flow  loss  coefficient  y*     is  converted  to  the  3*c- 
-ed  for  comparison  in  this  thesis  by 

%  =  h  7(2"dim) 

An  expression  for  }*  ^s  developed  that  is  a  function  of  the  circu- 
lation around  the  blades,  the  blade  width,  and  incompressible  flow, 

namely 

where  the  coefficient  K  depends  on  the  expansion  of  the  flow  in  the 
cascade.   The  value  of  K  is  determined  from  experimental  data  to  be 

Thus 


(17) 


In  Fig.    11   the   function    y<   is   shown   for  the   two-dimensional  efficiency 
of  Ref.    1, 

^         =  a>z   =  o.  9Z6Z 

/(2-dim)  'o 

Values  of  h/c  plotted  are  chosen  arbitrarily  but  Ref.  20  specifies 
that  the  derivation  holds  for  blades  that  are  not  excessively  short. 
There  is  assumed  to  be  some  two-dimensional  flow  in  the  channel  be- 
tween the  blades  at  the  center  section.   Experimental  data  are  shown 
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for  short  blades,  and  these  data  commence  at  [  — I  ■  2.0  and  go  to 
smaller  values.   So  the  plot  in  Fig.  11  is  based  on  eq.  (17)  and 
is  shown  for  a  value  of  I  — I  ■  2.0,  the  lower  limit  permitted.   For 
—  J  ■  1.33  and  1.25  the  mean  coefficients  are  about  107.  of  the  rep- 
resentative values  based  on  the  Vavra  equation  (3),  but  these  are 
specific  data  for  short  blades  not  arrived  at  using  the  equation 
above.   The  secondary  loss  as  computed  by  Markov  is  low  compared 
to  any  other  available  data.   When  complete  loss  coefficients  are 
computed  in  Ref .  20  the  values  are  comparable  to  any  other  method 
however,  since  the  tip  clearance  losses,  trailing  edge  losses, 
blade  height  corrections,  etc.,  used  bring  the  total  loss  to  the 
proper  level.   Any  use  of  the  Ref.  20  data  must  take  account  of 
this  fact. 

Of  interest  also  in  Ref.  20  is  a  relation  for  the  discharge 
flow  angle  of  a  turbine  blade  row.   The  relation  is 

A. 


(18) 


and  it  gives  more  accurate  results  than  the  other  widely-used 
equation 

.  A- 

A  (19) 

which  does  not  take  account  of  trailing  edge  thickness. 
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14.   The  tip  clearance  loss  as  represented  by  Ref.  21  Is  also 
ihown  in  Fig.  9.   The  equation  derived  is 

*" =  #!»+ * 

where  ££  -   .1011  +  4.667    (k/b) 

Tv  -  -.0422  +  2.790   (k/b) 

E 

h   ■  blade  height 
b   ■  blade  width 
k   ■  clearance 
and  the  results  are  shown  for  two  values  of  (h/b) . 

Also  given  in  Ref.  21  is  the  equation  for  the  minimum  induced 
drag  due  to  flow  in  a  clearance  gap  at  the  tip  of  an  airfoil  (as 
derived  by  Betz)  to  be 

3)~*z   &  f>  KFZ*  (21) 

where  K?  is  a  form  factor  depending  on  the  clearance,  blade  spacing 
and  relative  exit  angle  of  the  flow. 

Assuming  the  axial  component  of  velocity  to  remain  constant 
and  drawing  an  arbitrary  turbine  stage  velocity  triangle  (such  that 
all  angles  are  positive),  Fig.  12a,  the  circulation  (/"*)  can  be 
computed  by  integrating  around  the  closed  curve  in  Fig.  12b  to  give 


(5H*-  -"-) 


(22) 
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Then  the  lift  per  blade  is 

L  =f(-j)^L  -A  =  fxiAKfat  -M+< )       (23) 

The  lift  and  drag  forces  exerted  on  a  blade  can  be  seen  in  Fig.  13 
as  well  as  the  axial  and  peripheral  components  of  the  resultant 
force  F.   It  can  be  seen  that  for  each  blade  channel 

FA=Aj.(fi-fz)  (24) 

#  (25) 

The  secondary  flow  loss  will  be  associated  with  the  loss  from  the 
pressure  drop  in  the  axial  direction. 

To  arrive  at  a  loss  coefficient  expression  for  the  axial  force 
component,  it  is  necessary  to  define  the  loss  coefficient  in  terms 
of  the  pressure  loss,  as  shown  in  Fig.  14.   The  loss  coefficient  is 

defined  as 


*- 


7**4* 


Substituting 


%,-*%*  =  ffa-A) 
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and  for  incompressible  flow, 


*- =*+-£/" 


**     n  ^  2  jr  '  (27) 


and  collecting  terms, 


l*-*rir£-ir*-Z 


(28) 


Then  from  Fig.  13, 


/£  -  l-d^/f  +  P*+*., 


Using  the  above  derived  expressions  for  L  and  D  and  equating 


(29) 


(30) 


For  an  arbitrary  value  of  Kolhl  '  this  secondary  loss  is  plotted  in 
Fig.  15.   The  value  of  K?  \-r  )  can  be  adjusted  to  change  the  general 
inclination  of  the  family  of  curves  shown  to  any  desired  value.   That 
is,  a  lower  value  of  K-It-j  raises  the  curves  and  vice  versa. 
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15.  In  Ref .  22  an  attempt  is  made  to  find  a  mathematical 
solution  to  the  problem  of  secondary  flows  in  particular.  The 
approximate  relations  derived  are  not  able  to  satisfy  the  boundary 
conditions.  One  conclusion  drawn  from  the  derivations  is  that  the 
secondary  flows  are  a  direct  consequence  of  a  non-uniform  approach 
velocity,  and,  if  the  lift  is  assumed  proportional  to  the  square 
of  the  local  velocity,  then  the  secondary  flow  is  independent  of 
aspect  ratio.   These  conclusions  are  counter  to  the  ideas  of  most 
Investigators  in  the  field. 

The  tests  and  measurements  made  showed  odd  discrepancies  and 
it  was  felt  that  rotor-stator  induced  effects  may  have  been  the 
cause.   No  relation  was  deduced  which  is  directly  useful  here. 

16.  The  ideas  and  relations  presented  in  Ref.  23  and  24  were 
carried  through  and  refined  and  elaborated  in  Ref.  25.  All  three 
references  have  the  same  authors.  Therefore,  Ref.  25  will  be  dis- 
cussed here  as  being  most  useful,  and  parts  of  Ref.  23  and  24  will 
will  be  extracted  as  necessary. 

By  far  the  most  detailed  method  presently  available  for  deter- 
mining loss  coefficients  is  presented  in  Ref.  25.  The  loss  is 
broken  into  component  parts  similar  to  those  assumed  at  the  begining 
of  this  discussion  and  the  quantitative  data  offered  are  based  on 
tests  of  turbines. 
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The  loss  coefficients  of  Ref.  25  are 

P   -  P 
Rl    R2       Loss  of  total  head  pressure 


p       -   p  Total  pressure  at  blade  outlet-static  pressure  at  blade  outlet 

R^        A 


(32) 


The  loss  coefficients  used  in  this  thesis  are 


f  -/  -y  = 


'*j  —/ez 


and  the  necessary  conversion  is 


(33) 


<- 


(34) 


For  profile  losses,  two  figures  are  presented,  one  for  nozzle 
blades  and  one  for  impulse  blades.   These  are  reproduced  here  as 
Fig.  16a  and  Fig.  16b.   An  interpolation  method  is  proposed  for  rotors 
with  reaction  and  a  method  of  extending  the  profile  loss  to  conditions 
of  other  than  zero  incidence  is  also  shown. 

A  number  of  relations  for  the  loss  due  to  tip  clearance  are 

presented  in  Ref.  25  as 


r-"(i) 


(Stodola,  Reaction  turbine)    (35) 


(Meldahl)  (36) 


507.  Reaction  turbine)        (37) 


These   relations  are  shown   in  Fig.    9. 
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Expressions  are  also  derived  for  tip  clearance  loss  based  on  C 
and  c  if  the  ratio  spacing/blade  height  is  small  and  the  gas  turning 

Li 

angle  is  small.   For  our  uses  here,  the  gas  turning  angles  are  large. 
The  loss  relation  arrived  at  is 


(38) 


where  B  *  0.5  for  blades  with  clearance  and  B  ■  0.25  for  blades 
with  shrouds. 

Equation  (38)  is  used  to  find  the  range  of  loss  coefficients 
for  the  impulse  conditions  of  cos  0^   ■  1.00  and  cf.    ■  C(~.      Shown 
in  Fig.  9  is  the  band  of  values  corresponding  to  o( '     ■  0(7   ■  60  to  75  . 

The  data  available  to  Ainley  and  Mathieson  led  them  to  the  con- 
clusion that  the  effect  of  area  ratio  is  more  important  to  secondary 
flows  than  is  the  magnitude  of  the  gas  deflection.   The  area  ratio 

used  is 

A?     (annulus  area  at  outlet  from  blade  row)  (cos  0( '  ) 

— —     m        I,  a       (39) 

A.      (annulus  area  at  inlet  to  blade  row)  (cos  0    )  v   ' 

where  o(     »  gas  outlet  angle  measured  from  the  axial  and  3.   "  blade 
inlet  angle  measured  from  the  axial.  The  hub  ratio  (inner  diameter/ 
outer  diameter)  was  also  felt  to  exert  a  strong  influence.   This  led 
to  the  plot  of  (A2/A.)  /Cl  +  (ID.OD)]  versus  a  parameter  ^  used  in 
defining  the  drag  coefficient  due  to  the  secondary  flow  as 

cDs-KCLZl(*lc) 
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This  plot  is  reproduced  here  as  Fig.  17.   Ultimately  then,  a  loss 
coefficient  for  secondary  flows  is  presented  as 


K  -*\ 


<£*&  4?  A<*zo£, 


^i+i,o(j     &Ui*o(~ 


(41) 


Converting  these  values  of  /   to  ^  ,  a  plot  is  presented  as  Fig.  18 
for  a  wide  span  of  inlet  and  outlet  angles. 


For  trailing  edge  thickness  losses,  the  assumption  is  made  that 
the  previously  derived  losses  are  for  a  blade  of  trailing  edge  thick- 
ness which  is  27.  of  the  blade  pitch.   For  other  trailing  edge  thick- 
nesses a  correction  is  applied  according  to  a  plot  which  is  reproduced 

here  as  Fig.  19. 

G.   Results  of  Loss  Investigation 

The  results  of  the  loss  investigation  will  be  discussed  according 
to  the  original  breakdown  of  loss  components  mentioned  at  the  beginning 
of  this  section  of  this  thesis.   The  specific  turbine  design  used  in 
the  last  section  of  this  report  will  be  used  as  the  object  of  a 
quantitative  comparison  of  the  loss  relations  shown  here. 

1.   The  profile  loss  estimation  method  presented  by  Ainley  in 
Ref.  25  is  straightforward  and  complete  and  should  give  good  data 
for  the  "conventional"  blade  shapes  from  which  the  data  is  derived. 
An  anomaly  presents  itself  in  Fig.  10  however,  when  a  comparison  is 
made  between  the  profile  loss  coefficient  data  of  Ainley  and  the  data 
of  Markov.   The  blade  shapes  from  which  the  Markov  data  are  derived 
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are  very  similar  to  those  of  Fig.  25.   It  would  appear  that  such  blade 
shapes  would  have  higher  profile  losses  than  those  based  on  airfoil 
designs  which  usually  are  used  for  lower  flow  deflections.   The  data 
comparison  shows  this  is  not  the  case. 

Of  special  interest  in  the  Markov  data  in  Fig.  10  is  the  fact 
that  an  optimum  solidity  of  (s/c)  *  .625  to  .650  is  evident.   Higher 
solidities  apparently  block  the  flow  thereby  increasing  the  profile 
loss,  and  lower  solidities  encounter  separation  effects  since  the 
blades  act  more  like  individual  airfoils.   The  low  point  in  the  curve 
corresponds  to  guiding  the  flow  in  a  channel,  retarding  separation 
onset,  but  not  blocking  the  flow. 

2.   The  tip  clearance  losses  shown  in  Fig.  9  seem  to  give  a 
large  range  of  values  depending  on  which  curve  is  chosen.   The  sit- 
uation is  perhaps  more  reasonable  when  the  ratio  of  radial  clearances 
to  blade  height  is  restricted  to  realistic  values  less  than  .03. 
The  loss  estimate  from  the  Ainley  equation  is  obviously  outside  the 
flow  deflection  range  intended  for  that  relation.   The  other  curves 
for  tip  clearance  loss  for  reaction  blading  could  be  based  on  dif- 
ferent degrees  of  reaction  corresponding  to  the  fan  aspect  of  the 
curves,  i.e.,  high  reaction  blading  would  use  the  highest  curve. 
Some  reasonable  estimate  could  be  made.   The  only  relation  for 
impulse  blading  available  is  that  of  Markov.   The  Meldahl  relation 
is  the  only  one  which  accounts  for  the  trailing  vortex  inducing  a 
loss  even  though  the  tip  clearance  is  zero,  but  the  dependence  on 
blade  aspect  ratio  appears  unreasonably  strong  for  this  data  as  a 
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whole.   For  (h/b)  ■  1.00,  the  loss  at  (k/h)  ■  .02  is  greater  than  the 
highest  ordinate  shown. 

3.  The  Markov  data  is  all  that  is  available  for  making  an 
estimate  of  trailing  edge  thickness  loss  directly.   The  plots  in 
Fig.  20  indicate  the  trends  which  would  be  expected. 

The  Ainley  method  of  modifying  the  complete  loss  coefficient 
depending  on  the  relative  trailing  edge  thickness  seems  an  unfor- 
tunate way  to  present  the  data.  The  system  is  workable  but  indirect. 
It  would  be  more  useful  for  comparison  and  design  purposes  to  have 
the  trailing  edge  thickness  data  presented  directly. 

4.  The  situation  regarding  secondary  flow  loss  coefficients  is 
the  most  confused  of  all.   The  general  equation  by  Vavra  is  not  very 
precise  but  certainly  indicates  the  proper  trend  if  the  model  pre- 
sented in  Fig.  7  is  correct.   None  of  the  other  presently  available 
methods  investigated  will  even  show  the  proper  trends  for  high  flow 
deflections. 

The  Ainley  relation  plotted  in  Fig.  18  seems  to  have  depended 
too  heavily  on  stator  data  with  an  axial  entering  angle  and  low  flow 
deflections  in  the  rotor.   Otherwise,  the  curves  indicate  that  for 
a  given  entering  angle,  increasing  the  flow  deflection  decreases  the 
loss. 

In  Fig.  15  the  relation  derived  from  Meldahl  and  Betz  indicates 
the  proper  trend  for  low  flow  deflections  only.   Taking  the  plot  for 
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p  *  -45  as  an  example,  it  can  be  seen  that  for  flow  deflections 
greater  than  A0  a  90  the  indication  is  that  increasing  the  flow 
deflection  decreases  the  loss. 


The  Markov  relation  plotted  in  Fig.  11  suffers  the  same  deficiency 
pointed  out  for  the  Ainley  and  Meldahl  curves  plus  abnormally  low 
secondary  loss  coefficients  throughout  for  reasonable  blade  aspect 
ratios.   The  data  is  of  no  value  unless  used  specifically  with  Markov 
relations  for  all  losses.   In  Ref.  20,  when  sample  loss  estimations 
are  carried  out,  the  total  loss  coefficients  compare  well  with  the 
measured  values.   Apparently  the  profile,  tip  clearance  and  trailing 
edge  loss  coefficients  with  the  blade  height  corrections  used  in  the 
complete  method  counter  the  low  secondary  loss  estimates. 

5.   The  relations  now  available  can  be  used  to  make  a  quantitative 
comparison  of  the  component  losses.   The  turbine  designed  in  the 
following  section  will  be  used  as  a  model  to  show  the  comparisons. 
The  physical  dimensions  of  the  stator  and  rotor  are  as  listed  in 
Table  IV  and  these  are  used  in  the  relations  and  plots  so  far  developed. 
The  numerical  results  are  shown  in  Table  V  for  the  stator  and  Table  VI 
for  the  rotor. 
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Table  IV 
TURBINE  PHYSICAL  DIMENSIONS 


Stator 

Rotor 

Entering  angle 

0° 

69 . 60° 

Exiting  angle 

75° 

-71.57° 

Blade  height,  h 

1.165  in 

1.465  in 

Trailing  edge  thickness,  t 

.036  in 

.025  in 

Blade  chord,  c 

3.16  in 

1.355  in 

Blade  width,  b 

1.60  in 

1.375  in 

Spacing,  s 

1.80  in 

.9415  in 

Blade  thickness,  t 

.69   in 

.6485  in 

Tip  clearance  assumed,  (k/h) 

.00 

.02 

Table  V 
STATOR  LOSS  COEFFICIENT  COMPARISON 


Tip 

Trailing 

Author 

Profile 

Clearance 

Edge 

Secondary 

Loss 

Loss 

Loss 

Loss 

/avra,  Ref.l 

.070 

Ainley,  Ref.23 

.067 

Ainley,  Ref.25 

.047 

.041 

Markov,  Ref.20 

.030 

.032 

<.005 

Meldahl,  Ref.21 

.106 

.022 

Stodola,  Ref.3 

.040 
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Table  VI 
ROTOR  LOSS  COEFFICIENT  COMPARISON 


Author 

Profile 
Loss 

Tip 
Clearance 
Loss 

Trailing 
Edge 
Loss 

Secondary 
Loss 

Vavra,   Ref.l 

.231 

Ainley,  Ref.23 

~  .12 

Ainley,  Ref.25 

.148 

.13 

.168 

Markov,   Ref.20 

.068 

.025 

.032 

.025 

Meldahl,   Ref.21 

* 
.223 

.2601 

Stodola,    Ref.3 

.05* 

reaction 


As  can  be  seen,  the  component  losses  presented  in  Tables  V  and  VI 
can  be  combined  selectively  to  obtain  almost  any  desired  answer. 

The  loss  investigation  is  disappointing  in  that  no  clear-cut 
answers  are  provided.   The  complete  Ainley  method,  or  the  complete 
Markov  method,  may  produce  a  reasonable  estimate  but  the  correlation 
is  poor  at  best  for  component  losses. 

Any  desire  to  find  a  way  of  basing  a  new  design  on  component 
loss  considerations  is  presently  thwarted.   A  need  is  evident  for 
an  orderly  testing  program  to  find  the  magnitude  and  interaction 
of  component  losses. 
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IV.   Design  of  Turbine 

As  a  further  investigation  into  the  characteristics  of  a  turbine 
being  used  for  maximum  work  output,  a  single  stage  turbine  is  designed 
for  a  slightly  supersonic  velocity  leaving  the  stator  and  a  relative 
Mach  number  of  0.8  as  seen  by  the  rotor.   Then  the  overall  pressure 
ratio  Is  increased  and  the  off-design  performance  evaluated.   The 
number  of  calculations  required  are  prohibitive  for  hand  calculations 
so  the  equations  necessary  for  the  performance  evaluation  are  also 
programmed  in  Fortran  language  for  the  CDC  1604  Digital  Computer. 

A.   Loss  Coefficient  Data 

A  set  of  loss  coefficients  for  the  design  were  provided  by 
Vavra  from  some  unpublished  data.   These  loss  coefficients  are  shown 
in  Figs.  21,  22  and  23  and  will  be  used  for  the  design  and  the  per- 
formance e va lua  t  ions . 

The  abscissa  of  the  plot  of  stator  loss  coefficients,  M,  ,  is 

xs 

the  Mach  number  for  an  isentropic  expansion  from  the  total  pressure 
at  the  inlet  to  the  static  pressure  at  the  discharge  of  the  blade 
row.   This  can  be  seen  in  the  T-s  diagram  of  the  expansion  through 
a  turbine  shown  in  Fig.  24  to  be  the  isentropic  velocity  V  .  divided 
by  the  acoustic  velocity  based  on  T.  ,  .   When  reduced,  the  result  is 


'u**  /aw 


(42) 


-36- 


To  properly  determine  the  flow  area  through  the  stator,  the 
loss  coefficient  In  Fig.  20  noted  "for  flow  rates"  and  shown  in 
Fig.  24  as   3    is  used.   This  loss  coefficient  is  intended  to 
represent  the  losses  encountered  by  the  flow  in  the  channel  between 
blades,  but  is  not  intended  to  account  for  the  mixing  losses  and 
separation  losses  encountered  by  the  flow  leaving  the  blade  row. 
These  additional  losses  are  included  in  the  higher  loss  coefficients 
which  are  applied  to  velocity  determinations.   These  loss  coef- 
ficients are  noted  "for  velocities"  in  Fig.  21  and  as  f.    in  Fig.  24, 
Similar  considerations  are  used  for  the  loss  coefficients  through  the 
rotor.  ' 

The  loss  coefficients  are  corrected  for  blade  height  according 
to  the  curves  given  in  Fig.  23. 

B.   Flow  function 

For  steady  adiabatic  flow  the  stagnation  enthalpy  remains  con- 
stant along  a  streamline.   For  stators,  the  stagnation  enthalpy  is 


*  (43) 

and  in  Fig.  24,  to  a  different  scale,  H  would  be  represented  by 
T  and  h  would  be  represented  by  T,.   A  similar  situation  pertains 
to  a  relative  stagnation  enthalpy  defined  as 

Ht-^  +  g-Jt  +  kg 

for  the  condition  such  that  the  peripheral  velocity  at  the  entrance 
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of  the  rotor  equal  to  the  peripheral  velocity  at  the  exit  of  the 
rotor.  As  shown  in  Appendix  I,  the  mass  flow  rate  for  a  given  set 
of  inlet  conditions  can  be  expressed  by 


ait 


p.i      \p0i 
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where  n  is  the  poly  tropic  exponent,  P  is  the  total  pressure  ahead 
of  the  stator  blade  row  and  p  is  the  static  pressure  behind  the 
blade  row. 


As   shown   in  Fig.    24,   a   relative   total   temperature  and  relative 
total  pressure  can  be  defined  as 
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The  mass  flow  rate  through  the  rotor  is  then 
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Rearranging  the  above  expressions  for  flow  rate  in  non-dimensional 
terms  defines  a  flow  function, 
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and 
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It  will  be  necessary  to  investigate  whether  this  flow  function  can 
be  satisfied  for  each  set  of  conditions  imposed  across  a  blade  row. 
The  choked  flow  condition  corresponds  to  the  maximum  value  of  P/p 
which  the  blade  row  can  accomodate.   The  development  of  the  critical 
values  of  the  flow  function  and  pressure  ratio  can  also  be  found  in 
Appendix  I  to  be 
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A  computer  program  to  do  these  computations  is  explained  in  section  E 
below. 

C.   Geometry  and  Velocity  Triangles 

The  final  design  of  the  stator  and  rotor  blade  rows  is  as  shown 
in  Fig.  25  for  the  mean  diameter.   The  velocity  triangles  of  the 
average  flow  velocities  at  each  station  are  as  shown  in  Fig.  26. 

The  sign  convention  chosen  is  that  positive  angles  and  positive 
peripheral  velocity  components  point  in  the  same  direction  as  U,  the 
rotor  peripheral  velocity. 
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It  is  convenient  to  define  an  "exit  plane"  and  a  "discharge 
plane"  leaving  each  blade  row  as  shown  in  Fig.  25  by  the  (e)  and  (d) 
designations.   The  design  of  the  stator  and  rotor  is  such  that  the 
minimum  cross-sectional  area  of  the  channel  between  blades  occurs 
at  the  exit  plane  fe)  .   For  maximum  work  output,  the  pressure  drop 
across  the  stator  is  great  enough  to  cause  choked  flow  at  the  exit 
plane  plus  an  additional  expansion  of  the  flow  from  the  exit  plane 
to  the  discharge  plane.   The  conditions  for  such  after-expansions 
are  as  shown  in  Appendix  II  based  on  the  treatment  in  Ref.  26. 

In  Ref.  1  a  similar  development  is  shown  for  a  supersonic  flow 
entering  a  blade  row.   Since  the  flow  leaving  the  stator  is  super- 
sonic, it  is  possible  under  certain  conditions  for  the  rotor  to  see 
a  supersonic  flow  also  and  proper  account  must  be  taken  of  the  shock 
losses  involved. 

D.   Dimensionless  Parameters 

As  has  already  been  mentioned,  the  number  of  variables  to  be 
handled  is  greatly  reduced  if  the  mass  flow  expression  is  non-dimen- 
sionalized  to  form  a  flow  function  depending  only  on  pressure  ratio, 
specific  heat  ratio  and  the  polytropic  exponent  of  the  expansion 
process. 

In  Ref.  27  can  be  found  a  complete  development  of  dimensionless 
parameters  and  referred  values  for  minimizing  the  number  of  variables 
necessarily  handled  in  presenting  the  data  for  a  turbine  performance 
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analy8is.  Those  used  in  the  present  analysis  are  referred  rpm,  over- 
all pressure  ratio,  referred  flow  rate,  and  a  power  coefficient.  The 
necessary  coefficient  forms  are  computed  from 

Referred  RPM  -  RPM/  /~T~"  (51) 

Overall  pressure  ratio  ■  P  /p  (52) 

Referred  flow  rate  -  w  -/t""/P  (53) 

'  o  o 

Power  coefficient  -  HP/P  7  T  (54) 

of       o 

The  referred  flow  rate  and  power  coefficient  are  constant  in  this 

evaluation.      It   is  useful   to  refer  all   the  conditions  at  any  point 

in  the  turbine  to  inlet  conditions.   For  this  reason,  velocities  are 

carried  as  V/  V  T   and  W/  V  T   and  pressures  as  p/P  and  P„/P  . 
/o         /  o      r  ro      Ro 

E.   Design  Computations 

In  order  to  design  the  stator  for  choked  flow  at  the  exit  plane, 
the  flow  function  must  be  known  for  the  imposed  conditions.   A  com- 
puter program  was  written  to  compute  the  flow  function  for  a  given 
specific  heat  ratio,  2*  ,  pressure  ratio,  P/p,  and  loss  coefficient,  j  . 
As  a  by-product  the  program  also  computes  the  polytropic  loss  coef- 
ficient, X^   .   This  information  is  presented  in  Appendix  V  as  Table  III, 
Flow  Function  and  Polytropic  Loss  Coefficients,  for  X  from  1.25  to 
1.40,  ^  from  0.0  to  0.25  and  P/p  from  1.02  to  6.00.   The  particular 
flow  function  data  to  be  applied  here  for  ^  ■  1.37  are  also  included 
as  Fig.  27.   The  computer  program  is  included  in  Appendix  III. 
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The  following  data  were  assumed  given  for  the  design: 


Mass  flow  rate 

■ 

100.  lbm/sec 

Total  pressure  at  inlet 

- 

1000  psia 

Total  temperature  at  inlet 

s 

1460°R 

Specific  heat  ratio 

- 

1.37 

Gas  constant,  R 

- 

421.5  ft  lbf 

Mean  diameter  of  stator 

m 

23.0  in 

Mean  diameter  of  rotor 

m 

23.0  in 

Design  rpm 

- 

13000 

The  stator  blade  design  was  provided  by  Vavra  from  unpublished 
data  and  is  as  shown  in  Fig.  25.   It  is  now  necessary  to  select  the 
blade  height  and  number  of  blades  so  the  choked  flow  condition  at 
the  exit  is  satisfied. 

A  blade  height  of  1.16  inches  is  taken  as  a  first  estimate.   Then, 

from  Fig.  21  the  loss  coefficient  for  flow  areas  at  M.   ■  1.00  is 

is 

found  to  be   5    *  .053.   From  Fig.  23  the  blade  height  correction 

'  aO 

is  K,  ■  0.872  and  so 


f-. 


872  (.053)  -  .0462 


Using  linear  interpolation  on  Fig.  27  for  this  loss  coefficient 
between  the  maximum  values  of  (p   ,  i.e.,  (p      ,  it  is  found  that 

0  _   -  .65811 


cr 


(!)    -».« 


cr 
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Recalling  the  flow  function  equation  (48),  it  is  found  that 


0        -  ^dK,f^     *   13 

A  =       13.82  2 

7*>ea4r  .65811  u  ln 


From  Fig.  25  the  measured  minimum  distance  between  blades  at 
the  exit  plane  is  a  »  .448  in.  and  s  ■  1.8  in.  is  given.   The  number 
of  blades  is  then 


/ 


yJL;    -  40.2 

*  (55) 


In  order  to  prevent  harmonics  in  the  wake  pattern  from  the  stator 
as  seen  by  the  rotor,  a  prime  number  of  stator  blades  is  chosen. 
In  this  case,  z  ■  41  blades  will  suffice.   For  41  blades, 


4   =r  ^"(23.0)   -  1.763 

41 


so  all  measured  values  from  Fig.  25  must  be  scaled  by  the  scale 

factor, 

1    763 
Scale  Factor  »        •    g     -  0.9794 

The  blade  height  must  now  be 

K,    "    —     *  TT-TrTTTfl Q7Q/x      ■  1.165   in  (56) 

SI    za    41.0  (.448  x  .9794)  v 

Recomputing  (b        with  a  new  loss  coefficient  based  on  this  new  blade 
height  makes  no  measurable  difference  in  exit  area  required  so  the 
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stator  row  design  is  now  fixed.   To  simplify  the  iteration  on  (J  a 
computer  program  was  written  to  find  the  choked  flow  values  of  <fi      , 
(P/p)   and  ^5  .   This  computer  program  is  included  in  Appendix  III 
and  the  computed  data  are  presented  as  Table  VII  in  Appendix  V. 

The  flow  conditions  at  the  stator  exit  plane  are  determined 
from  Fig.  24  as  follows: 

-f*    «-     1       -   .538 

p  1.86 


8458 


1470 


JL-JiZ£±jK  -121.52 
Z      I    *-'      % 

When  a  check  is  now  made  of  the  Mach  number  at  the  exit  plane  in 
order  to  check  for  choked  flow  it  is  found  that 

\    =    TJ±=        "  0.97 
fit**. 

This  computed  Mach  number  is  the  average  Mach  number  of  the  flow  at 

the  minimum  area.   The   highest  Mach  number  at  this  area  is  equal 
to  1.00  but  it  is  averaged  with  the  Mach  number  of  the  flow  in  the 
boundary  layer  as  accounted  for  by  the  loss  coefficient.   In 
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Appendix  IV  the  derivation  of  the  Mach  number  for  choked  flow  in  a 
polytropic  expansion  is  presented.   The  "critical  average"  Mach 
number  is  found  to  be 


"*      '/+&(*-')       I  *-'  (IV' 9) 


This  condition  will  also  be  encountered  later  in  the  computer  program 
for  the  flow  through  the  rotor  and  proper  account  must  be  taken  of 
it  there. 


Using  the  relation  from  Ref.  20  for  flow  exit  angle  and  the 

Bl 


measured  blade  exit  angle  from  Fig.  25  of  0(n,    •  73.9  gives 


On*,    -  ^ =      .268 

(18) 


*i   «  75° 


For  the  rotor,  the  blade  height  is  permitted  to  be  approximately 
0.30  inches  greater  than  the  stator  blade  height  in  order  to  permit 
the  streamlines  to  expand  and  offer  no  resistance  to  flow  leaving  the 
stator  but  not  have  dead  areas  above  and  below  the  useful  flow  area 
of  the  rotor  blade.   The  rotor  blade  height  will  be  taken  then  as 
h   -  1.465  in. 

In  order  to  satisfy  the  requirement  for  the  rotor  to  see  a 
relative  Mach  number  of  0.8,  some  after-expansion  is  required.   In 
order  to  evaluate  the  proper  amount  of  after-expansion,  the  equations 
in  Appendix  II  were  coded  for  the  computer  to  solve  for  any  applied 
pressure  ratio  between  the  exit  and  discharge  planes  for  a  given  "3^  , 
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Mach  number  and  blade  angle.  This  computer  program  Is  included  in 

Appendix  III  and  the  solution  for  1f  »  1.37  and  o( .    ■  75  over  a 

range  of  pressure  ratios  is  presented  here  as  Pig.  28.   The  computer 

program  can  be  used  for  any  specific  hsat  ratio,  blade  angle,  Mach 

I 

number  and  pressure  ratio  range. 


t  t 


After  some  iteration,   it  is  determined  that  p./P     ■  0.460  should 

1  o 

give  approximately  the  correct  pressure  ratio  across  the  expansion 
area  from  the  exit  plane  to  the  discharge  plane, 

!l  .  PJL    .    Jgl    .  o.„5 
Pe     Pe     .538 

Now  from  Fig.    28  it   is   found  that    &Q(  *  0.2°     and     *J      *  .0008.     The 
total   loss  coefficient  to  be  appliss'   is  based  on  a  blade  height  of 
1.165  inches   for  the  stator.     lfcars  is      ^  •   .1045  from  Fig.    21  and 
K,    ■   .838   from  Fig.   23  for  velocity  computations.     Therefore 

f,  -  fCKc)  +  £  -  °883 

Now  the  velocity  leaving   the  stator  is  determined  with  Fig.    24 

as   reference 

&Z      =    (Z_  _  T,^\(/  -  f,)    -  .1892  (.9117)  •  .1725 
\T.  Tt    J 


T0 


•ft         I     V-/       Te 


131.60 
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To  find  the  relative  angle  and  relative  velocity  of  the  flow,  it  is 
necessary  to  construct  the  velocity  triangle  as  shown  in  Pig.  26. 


it  7zo< 


•  34.18 


o 


«;    -    75°  -  0.2°  •  74.8 

vT  rT 

Vu±-1L4*~    74.8°.  l270 
*JM-      *&     -<L     -,2.82 

fT     frT     /T 

3=fe : — S"      ■  99.07 


60° 


As  a  check  on  the  requirement  of  a  relative  Mach  number  of  0.8  as 
seen  by  the  rotor 

M.   =    K"?   ..fLfll  .080 


To  design  the  rotor,  several  considerations  must  be  kept  in 
mind.  The  flow  deflection  should  be  as  great  as  is  practicable  for 
maximum  work  output.  The  trailing  edge  thickness  will  be  arbitrarily 
selected  as  .025  inches,  and  the  spacing  and  blade  design  will  be 
chosen  so  the  minimum  cross  section  area  of  the  flow  passage  between 
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• 


. 


•    . .' 


blades  will  occur  at  the  "exit"  plane.  The  rotor  blade  width  is 
selected  as  1.375  inches  and  impulse  conditions  only  will  be  used 
for  unchoked  flow,  i.e.,  p  ■  p.. 

Recalling  that  a  relative  total  temperature  and  relative  total 
pressure  are  defined  as 


-IB.  m    — .  y. '       „    m  —  +  — iJ £--.9253 

&  _  *  /  5/  X&     „„ 

«    — --/ /  •  .6955 


The  flow  function  as  applied  to  the  rotor  is 


(46) 


(47) 


^  «  ^cfZJZ 


if 


Az    °g.i  9    f  (48a) 


and  this  expression  will  be  used  to  determine  the  exit  area  of  the 
rotor.   For  the  prescribed  impulse  conditions 

1*1  -  -6955    .  S1<. 
PJ?  T4600  "  l'515 

and  the  value  to  be  used  on  the  abscissa  of  the  plots  for  rotor 
loss  coefficients  in  Pigs.  22  and  23 


-A       ■  SSA  ^°       ■       2SQS 

Vx  131.61  *"98 
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.... 

■  ■  ■  .-.    ■     .   ■■ 


With  this  value  the  rotor  loss  coefficient  for  flows  Is  found  to  be 

V j  •  .819  (.192)  •  .1572 

Now  using  Fig.   27  for  this  pressure  ratio  and  loss  coefficient  and 
using  linear  interpolation  the  flow  function  is 


P   -  .58953 
and  from  eq.    (48a)  above 

^A2     •  19.13 


A2     •  32,44     in2 

An  impulse  rotor  operating  under  conditions  as  these  should  have 

a   ratio 

It)"*' 

and   so 

s  •  0.7  (1.375)  -  0.962 

For  the  mean  diameter  of  23.0  lnchea,  the  number  of  blades  on  the 
rotor  is  then 

Z2  "      '^1 962°  ^  *  75,1  blades 

Since  an  odd  number  of  blades  was  chosen  for  the  stator  the  rotor 
should  have  an  even  number  of  blades  to  reduce  the  possibility  of 
resonance  induced  vibrations  and  the  resulting  metal  fatigue. 
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Arbitrarily  choosing  76  blades 


.  .  ■%("■<>)  .  .9415  m 


and   Che   distance   between  biadee,   a,    to  obtain   the   required  area  of 
32.44  square  Inches  Is 

A2 
a  -    „     £        -  .273  in 

The  exiting  blade  angle  can  be  seen  from  Fig.  25  to  be 


i4 


A- 


•  .3162 


/«  "  71S7< 


The  blade  design  Itself  Is  based  on  a  circular  arc  contour  of 
the  concave  face.  The  geometric  relations  of  the  design  have  been 
programmed  by  Vavra  for  computer  solution.   Using  this  computer 
program  and  requiring  solutions  for  74,  76,  78,  80  and  82  blades, 
the  printout  of  the  computer  solution  Is  presented  as  Fig.  29.  The 
values  necessary  to  draw  the  blade  design  are  shown  in  Fig*  30  and 
have  been  used  to  draw  Fig.  25. 

The  solution  for  76  blades  appears  valid  and  so  the  other  solutions 
are  not  used.  The  resulting  It-}  ■  0.691.  A  check  on  this  spacing  can 
be  made  by  Brillngs  Rule,  an  early  rule-of- thumb  for  steam  turbines 
which  is  found  in  Re f.  3  to  be 
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( j)  •  2.5  sin  Ifi  *   where  A*   -  90  -  A^  (57) 

j-|j  •  2.5  sin  36.86°  •  1.5 

(f)-!1?  -6667 

Also  Ref.  28  presents  a  relation  for  spacing  and  blade  angles 
based  on  an  aerodynamic  load  coefficient  lr_  which  is  found  to 
have  a  nearly  constant  value  over  a  wide  range  of  turbine  designs. 
The  relation  is 

and  so,  using  proper  angles  as  defined  in  Ref.  28  for  this  case 


(58) 


40° 

,83° 


/si  m     0.8  sin  20.40° 

I*'   "     20  sin  18.43°  sin  38. 1 


The  relations  by  Brillng  and  Zweifel  show  at  least   that   the  choice 
of  76  blades   is  reasonable  for  this  design.     The  relation  of  Zweifel 
in  particular  is  based  on  two-dimensional  considerations  and  does  not 
account   for  blade  height.      Some  unpublished  test  data  available  to 
Vavra   indicates   the  value  of  0.691  is  a  good  one  for  this   size   turbine. 

To  complete   the  computations  necessary  to  determine  the  power 
output  of  the   turbine   for  a  given  set  of  conditions,    the  velocity 
triangle  after  the  rotor  is  computed. 
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For  no  after-expansion,  y^2     »    y^po  ■  -71.57   ,   and  finding  loss 


B2 
coefficients  in  the  manner  previously  described, 


?Z<U  _    7g//^    \  *  .   .9263  I  - 


'    460°  I"2701  -.8275 


6955 


^=/A'_W>-&)..ono 


-2-  *     -21 -£.    •  .8543 


-M.   »  JlU*-*~    .34.3, 

7T    /  «•-/    £ 


Using  relations  as  before  with  Fig.  26  as  reference,  it  is  found 

aS4      =*  -  45.93 
and  the  power  output,  with  the  proper  conversion  factors   is 


HP  "     gj     (U1  Vul  "  U2Vu2)   (1'414>  "  48»682     HP  <59> 


These  computations  are  verified  by  the  computer  program.      The   first 
page  of  the  printout  of  the  computer  program  is  presented  as  Fig.    31. 
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V'  Computer  Program 

A.  General  Considerations 

The  computer  program  for  the  evaluation  of  the  turbine  performance 
is  Intended  to  be  as  general  as  possible  within  the  framework  of 
the  turbine  design  method  used.  Any  single  stage  turbine  so  designed 
can  be  evaluated  by  entering  new  physical  dimensions,  entrance  con- 
ditions and  loss  coefficient  data.  The  static  pressures  to  be  im- 
posed ahead  of  and  behind  the  rotor  must  be  programmed  by  the  user 
within  the  main  program.  The  loss  data  must  be  of  the  same  format 
as  Pigs.  21,  22  and  23. 

B.  Main  Program 

1.   The  main  program  is  intended  to  make  the  major  computations 
and  comparisons  and  provide  control.   The  subroutines  are  provided 
to  do  the  lengthy  iterative  processes  and  repetitive  calculations. 
The  flow  charts  for  the  main  program  and  the  subroutines  are  in- 
cluded with  the  program  listings  found  in  Appendix  III. 

The  basic  concept  of  the  computer  program  is  as  follows: 

a.  With  the   input  data  given,   compute   the  design  condition 
including  the  after-expansion  at  stator  discharge. 

b.  Reduce  the  static  pressure  through  the  impulse  rotor  in 
increments,  computing  velocity  triangles  and  power  at  each  pressure 
increment.  The  shock  condition  for  supersonic  flow  entering  the 
rotor  is  computed.   Also  the  flow  function  through  the  rotor  is 
checked  for  the  choked  condition. 
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c.  When  the  choked  condition  through  the  rotor  is  attained, 
the  static  pressure  at  the  rotor  discharge  is  lowered  in  increments 
and  the  after-expansion  condition  computed.   The  resulting  power  and 
velocity  and  thermodynamic  conditions  are  computed  for  each  pressure 
increment. 

d.  The  process  is  terminated  when  the  loss  coefficient  due 
to  the  rotor  after-expansion  is  reduced  to  zero. 

The  program  will  also  terminate  if  the  loss  coefficient  due  to 
after-expansion  out  of  the  stator  goes  to  zero  or  if  the  loss  coef- 
ficient due  to  an  entry  shock  into  the  rotor  goes  to  zero.   Any  con- 
tinuation of  the  program  beyond  these  limiting  conditions  would  lead 
to  conditions  which  violate  the  Second  Law  of  Thermodynamics,  as 
explained  in  Ref.  1. 

The  loss  coefficient  reaches  a  maximum  and  then  declines  to  zero 
or  negative  values  as  the  imposed  pressure  ratio  across  the  control 
area  is  decreased  in  increments  for  the  after-expansion  case.   The 
exact  reason  for  this  is  not  known.   A  check  was  made  to  see  if  the 
axial  component  of  the  Mach  number  of  the  flow  leaving  the  blade 
row  reached  M  ■  1.00  when  the  loss  coefficient  reduced  to  zero.   It 
is  found  this  is  not  exactly  the  case  but  that  the  axial  component 
of  the  Mach  number  is  less  than  one.   Further  work  is  needed  on  this 
point. 

2.  Some  pertinent  comments  may  be  useful  about  the  main  decision 
points  in  the  main  program.  The  flow  chart  in  Appendix  III  shows  that 
the  first  check  is  for  the  magnitude  of  the  static  pressure  between 
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the  stator  and  rotor  compared  to  the  static  pressure  in  the  exit  plane 
of  the  stator.   For  a  stator  designed  as  a  converging-diverging  nozzle, 
the  Hach  number  at  the  exit  plane  may  be  greater  than  M   ■  1.00.   If 
the  static  pressure  between  the  stator  and  rotor  is  greater  than  the 
static  pressure  at  the  stator  exit  plane,  the  shock  condition  is 
computed  at  the  stator  exit.   If  the  pressures  are  equal,  there  is 
no  flow  deflection,  and  if  the  static  pressure  between  the  stator 
and  rotor  is  lower  than  the  exit  plane  pressure,  the  after-expansion 
case  is  computed. 

The  next  decision  is  based  on  the  relative  Mach  number  as  seen 
by  the  rotor.   For  the  supersonic  case,  the  shock  of  the  flow  in  the 
rotor  entry  is  computed.   For  subsonic  flow  the  component  of  relative 
velocity  in  the  direction  of  the  rotor  blade  entering  angle  is  used 
to  compute  flow  conditions. 

The  flow  function  and  critical  pressure  ratios  are  used  to 
determine  choked  conditions  in  the  rotor.   As  can  be  seen  in  Appendix  IV 
the  Mach  number  of  the  flow  at  the  rotor  exit  plane  is  an  average 
value,  and  the  critical  Mach  number  varies  with  the  loss  coefficient 
as  shown  there.   Since  the  impulse  conditions  for  the  rotor  are 
specified,  some  care  is  taken  by  an  iteration  process  to  find  the 
p.  •  p~  condition  for  choked  flow  before  the  program  continues. 

Once  the  choked  condition  in  the  rotor  is  attained,  the  con- 
ditions forward  of  the  rotor  exit  plane  remain  fixed.   The  Mach  number 
of  the  flow  at  the  rotor  exit  plane  is  set  at  M^   »  1.00  and  the 
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static  pressure  behind  the  rotor  is  lowered  in  increments.   The  after- 
expansion  behind  the  rotor  is  then  computed  until  the  limiting  con- 
dition is  reached. 

C.   Subroutines 

1.  The  subroutine  called  TRNGL  computes  the  velocity  triangle 
based  on  rpm,  absolute  velocity  and  direction  or  based  on  rpm, 
relative  velocity  and  direction. 

2.  The  subroutine  called  AFTER  computes  the  deflection  of  the 
flow  due  to  shock  or  expansion  of  the  flow  out  of  a  blade  row.   The 
flow  angle  is  incremented  from  the  blade  angle  and  the  pressure  ratio 
is  computed  and  compared  to  the  imposed  value.   The  process  con- 
tinues until  the  flow  deflection  is  found  which  corresponds  to  the 
pressure  ratio  imposed.   The  detailed  equations  are  derived  in 
Appendix  II,  based  on  Ref.  26. 

3.  The  subroutine  for  computing  the  flow  function  for  a 
pressure  ratio  imposed  across  the  blade  row  is  named  CPHI.   The 
polytropic  loss  coefficient  corresponding  to  the  design  flow 

loss  coefficient  is  determined  by  iteration.   Then  the  flow  function 
to  satisfy  the  imposed  conditions,  critical  flow  function,  and 
critical  pressure  ratio  are  computed  from  the  equations  developed 
in  Appendix  I.   The  pressure  ratio  was  found  to  be  more  sensitive 
in  fewer  significant  figures  than  the  flow  function  itself  so 
decisions  in  the  main  program  are  based  on  pressure  ratio  comparisons 
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4.  Subroutine  BEFORE  is  used  to  compute  rotor  entry  shocks  based 
on  the  relations  in  Ref.  1  which  in  turn  are  based  on  considerations 
similar  to  those  of  Ref.  26  for  after-expansion.   Since  the  flow 
deflection  is  known,  an  iteration  process  is  not  required  and  the 
computations  are  straightforward.   For  cases  which  are  not  directly 
solvable  by  the  basic  equations,  the  component  of  the  velocity  in 

the  direction  of  the  blade  angle  is  taken  as  the  useful  resultant. 

5.  The  CURVE  subroutine  is  written  as  a  means  of  finding  a 
data  point  on  a  curve.   The  loss  coefficient  curves  "for  velocities" 
from  Figs.  21,  22  and  23  are  stored  in  the  computer  in  the  form  of 
arrays.   The  subroutine  writes  the  second-order  equation  of  the 
appropriate  three  stored  data  points  nearest  to  the  input  abscissa, 
and  then  finds  the  ordinate  on  the  curve  with  the  input  abscissa. 
The  subroutine  is  arranged  in  such  a  manner  that  the  stored  data 
must  have  equal  abscissa  increments  and  the  first  abscissa  point 
must  be  zero.   This  means  arbitrarily  extending  a  curve  such  as 

Fig.  21  to  M.   "0.0  whether  the  data  at  low  values  of  M.   will  be 

■         is  is 

used  or  not. 

D.   Results  of  Computer  Program 

1.   The  turbine  design  data  given  in  the  preceding  section  of 
this  thesis  were  entered  into  the  computer  program  and  the  resulting 
printout  for  the  design  condition  only  is  presented  as  Fig.  31.   The 
pressure  ratio  across  the  turbine  was  then  increased  in  increments 
(by  lowering  p.  and  p_  while  holding  P  constant)  and  the  resulting 
computations  are  included  as  Table  VIII  in  Appendix  V.   The  rpm  was 
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also  decreased  from  13000  to  9000  in  1000  rpm  increments  and  the 
complete  pressure  range  imposed  at  each  rpm.   These  results  are 
also  included  in  Table  VIII. 

2.  The  non-dimensional  parameters  previously  discussed  were 
then  plotted  in  order  to  make  cross  plots  which  ultimately  result 
in  a  portion  of  the  turbine  performance  map. 

3.  The  data  immediately  available  from  the  computer  printout 
permits  the  direct  plotting  of  Fig.  35,  Turbine  Efficiency  vs. 
Pressure  Ratio,  with  the  rpm  as  a  parameter.   The  anticipated  gradual 
loss  of  efficiency  with  increasing  pressure  ratio  is  evident.   Of 
interest  is  the  change  in  slope  of  the  curve  at  the  point  where 
rotor  after-expansion  commences.   The  indication  is  that  the  increase 
in  velocity  due  to  after-expansion  causes  the  peripheral  component 

of  the  velocity  to  increase,  even  though  a  deflection  of  the  flow 
toward  the  axial  direction  also  occurs.   Therefore,  the  rate  of  loss 
of  efficiency  is  decreased,  although  the  change  is  very  small.   A 
reduction  in  rpm  contributes  to  a  large  loss  in  efficiency  at  any 
given  pressure  ratio. 

4.  In  Figs.  36  and  37  the  most  important  indication  is  as 
follows:   Fig.  37  shows  that  for  a  fixed  design  at  a  pressure  ratio 
of  (P  /p»)  -  2.174,  the  increase  in  power  for  an  increase  in  overall 
pressure  ratio  is  direct  and  marked  up  to  a  pressure  ratio  of  about 
4.5  or  5.0.   Further  large  increases  in  pressure  ratio  net  very 
small  gains.   In  Fig.  36  the  corollary  is  shown  for  efficiency.   The 
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initial  increase  in  pressure  ratio  and  power  output  causes  a  certain 
loss  of  efficiency  but  the  remaining  increase  in  pressure  ratio  and 
power  output  causes  a  drastic  loss  of  efficiency. 

An  interesting  comparison  can  now  be  made  with  the  high  head 
coefficient  data  from  the  first  section  of  this  thesis.   In  that 
section,  each  point  on  each  curve  represents  a  different  turbine 
design.   The  conclusion  was  that  if  a  required  high  power  output 
transcends  the  need  for  high  efficiency,  then  a  very  high  head 
coefficient  should  be  used  for  the  turbine  design.   The  initial 
increase  of  head  coefficient  from  the  optimum  value  causes  a  large 
loss  of  efficiency.   Once  this  initial  efficiency  loss  is  accepted, 
further  increase  in  head  coefficient  require  little  more  loss  of 
efficiency.   Now  in  this  section  concerned  with  the  off-design  per- 
formance of  a  particular  fixed  turbine  design,  the  conclusion  is: 
additional  power  can  be  obtained  at  a  higher  pressure  ratio  than 
design  but  only  in  the  initial  increase  of  pressure  ratio.   The 
extremely  high  loss  in  efficiency  is  encountered  if  the  last  possible 
increment  of  power  is  demanded. 

5.   The  turbine  performance  map  presented  as  Fig.  38  serves  to 
indicate  the  pertinent  trends  for  the  parameters  investigated.   The 
lines  of  constant  efficiency  form  only  a  portion  of  the  elliptical 
curves  they  become  over  the  complete  range  of  referred  rpm  and 
power  coefficient.   The  only  conditions  investigated  in  this  thesis 
are  at  high  power  coefficients.   The  peak  efficiency  for  this  turbine 
would  occur  in  the  lower  right  hand  corner  of  Fig.  38,  i.e.,  in  a 
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lower  power  coefficient  range  than  shown  here  and  in  the  referred 
rpm  range  of  about  350. 

6.   This  turbine  designed  for  high  head  coefficients  can  now  be 
compared  to  the  theoretical  values  predicted  in  the  first  section  of 
this  thesis.   The  head  coefficient  is  defined  in  eq.  (1).   The 
isentropic  velocity  C  is  given  by  reference  to  the  ^T.   of  Fig.  24  as 


r 


zfscr 


(60) 


For  pressures  and  temperatures  at  the  design  point 


■4-  = 


c0  /jz 


-iZ 


137.8391 
34.1441 


16.2 


(1) 


Referring  now  to  Table  I  at  0{.    -  75°,  r  =0.0  and  k.,   =  16.0  the 

1  is 

values  determined  from  the  theoretical  conditions  can  be  compared 
to  the  turbine  design.   The  comparison  is  shown  in  Table  IX  in  the 
columns  Theoretical  #1  and  Turbine  Design. 


The  ratio  V  0/V  ,  was  arbitrarily  selected  equal  to  one  for  the 
m2  ml  n 

theoretical  development.   This  was  not  a  selected  parameter  for  the 

turbine  design  and  the  discrepancy  between  the  two  values  of  V  „/V  . 

m2  ml 

are  the  basis  for  the  differences  seen  in  Table  IX.   Since  the  ratio 

V  ,,/V  .  was  not  required  to  be  equal  to  one  in  the  turbine  design, 
mZ  ml 

the  flow  is  deflected  more  in  the  rotor.   This  results  in  higher 
absolute  values  for  u~   and  °»o-   Since  the  flow  deflection  is 
greater,  the  peripheral  component  of  the  absolute  velocity  (V  ~)  is 
greater  and  the  power  output  and  efficiency  are  increased.   Otherwise, 
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the  theoretical  predictions  are  good.   The  V  0/V  n  value  from  the 

m2  ml 

turbine  design  can  be  entered  into  the  equations  and  the  corres- 
pondence should  be  more  exact.   This  was  done  and  the  resulting  data 
is  the  column  headed  Theoretical  #2  in  Table  IX. 

Table  IX 
COMPARISON  OF  THEORETICAL  AND  TURBINE  DESIGN  VALUES 


*1 

Theoretical  #1 

Turbine  Design 

Theoretical  #2 

75° 

74.8   (after-expansion) 

75°    * 

h 

0     * 

0 

0     * 

R2/Rx 

1.0   * 

1.0 

1     * 

V  ,/V  - 
m2  ml 

1.0   * 
.95  * 

.784 
.955 

.784  * 
.95  * 

rR 

.95  * 

1.00 

.95  * 

f 

.849 

.852 

.833 

* 
r 

0.0   * 

0.0 

0.0   * 

kis 

16.    * 

16.2 

16.2   * 

7 

"e 

.60111 

.723 

.60358 

2.263 

2.42 

2.031 

vmi/ui 

.984 

.998 

.990 

W2/Wl 

.827 

.852 

.812 

i 

69.78° 

69 . 60° 

69.82° 

-65.30° 

-71.57° 

-70.54° 

«2 

-49.17° 

-59.84° 

-57.02° 

selected,  not  computed 
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7.  The  curves  of  Fig.  26  are  shown  again  in  Fig.  39  with  all 
values  referred  to  the  design  point.   All  values  are  then  shown  as 

a  percent  of  the  design  point  values.   Fig.  39  demonstrates  the  strong 
effect  of  rpm  on  power  output  and  efficiency.   At  707.  design  rpm,  it 
is  not  possible  to  develop  design  power  output.   At  about  907. 
design  rpm,  the  design  power  can  be  developed  at  957o  of  design 
efficiency.   A  maximum  power  of  about  1257.  of  the  design  value  is 
the  maximum  which  can  be  developed  at  907.  design  rpm.   At  the  design 
rpm,  1257.  design  power  can  be  developed  with  a  loss  of  efficiency 
of  107..   When  greater  power  increments  are  required,  the  efficiency 
loss  becomes  increasingly  large.   A  407.  increase  in  power  can  be 
developed  but  at  only  72.57.  of  design  efficiency. 

At  design  rpm,  1257.  design  power  requires  an  overall  pressure 
ratio  of  3.5,  and  1407.  design  power  requires  an  overall  pressure 
ratio  of  8.5. 

8.  As  this  thesis  was  being  published,  it  was  determined  that 
the  computer  program  finds  the  off-design  performance  of  a  turbine 
in  which  impulse  conditions  are  maintained  until  choked  flow  occurs 
in  the  rotor.   A  rotor  blade  height  of  1.648  inches  is  required  to 
satisfy  such  a  condition.   When  this  new  rotor  blade  height  is 
introduced  into  the  computer  program,  the  only  direct  influence  is 
on  the  blade  height  correction  factor  for  loss  coefficients.   The 
performance  computations  change  is  negligibly  small. 
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VI.   Conclusions  and  Recommendations 

A,   Conclusions 

1.  If  a  requirement  for  a  high  shaft  power  output  turbine  is 
encountered  for  which  weight,  size  and  space  limitations  preclude 
having  several  stages  operating  at  peak  efficiency,  then  a  slightly 
higher  head  coefficient  and  slightly  lower  efficiency  compromise 
design  should  not  be  made.   Designing  for  power  outputs  higher  than 
optimum  should  be  based  on  extreme  power  output  per  stage.   The 
initial  loss  of  efficiency  is  great  when  a  slightly  higher  than 
optimum  head  coefficient  is  considered.   Not  much  more  efficiency  is 
lost  by  imposing  extreme  values  of  head  coefficient. 

2.  From  the  available  literature,  no  straightforward  method 

of  computing  losses  is  found  which  correlates  well  with  the  available 
data.   This  appears  to  be  a  fruitful  area  for  future  research  and 
testing. 

a.   The  available  data  indicates  the  tip  clearance  loss  is 
not  zero  even  for  impulse  conditions  and  no  clearance  gap.   The 
trailing  vortices  induced  at  the  ends  of  the  blades  introduce 
velocity  components  which  do  no  useful  work.   Clearance  gaps  within 
the  annulus  boundary  layers  should  give  tip  clearance  loss  values 
about  the  same  as  zero  clearance.   As  the  clearance  gap  is  increased, 
the  loss  increases.   Also  the  loss  increases  as  the  degree  of  reaction 
is  increased.   A  high  degree  of  reaction  means  a  high  pressure  drop 
across  the  blade  row  and  such  a  pressure  differential  is  conducive 
to  high  leakage  across  the  blade  tips. 
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The  loss  due  to  tip  clearance  is  also  affected  by  the  relative 
size  of  the  gap  to  the  blade  height.   The  effect  of  leakage  across 
the  tips  of  the  blades  on  the  secondary  flow  will  make  measurements 
difficult  to  separate  into  tip  clearance  and  secondary  flow  losses. 

b.  The  profile  loss  is  a  function  of  blade  thickness  to 
chord  ratio,  spacing,  flow  deflection  and  probably  general  blade 
shape.   This  loss  offers  the  most  possibilities  for  simple  and 
accurate  measurements  in  a  cascade  test  rig. 

c.  The  trailing  edge  thickness  losses  are  probably  in- 
sensitive below  a  certain  ratio  of  trailing  edge  thickness  and 
spacing,  and  then  assume  a  real  importance  at  higher  values.   These 
losses,  too,  should  be  relatively  simple  to  measure. 

d.  The  secondary  loss  measurement  offers  the  most  difficulty. 
No  known  theory  accurately  predicts  the  loss.   Indications  are  that 
the  secondary  loss  is  primarily  dependent  on  the  flow  turning  angle 
and  ratio  of  blade  height  to  the  chord  or  axial  width  of  the  blade. 
The  tip  clearance  effect  is  also  of  primary  importance  and  will  be 
difficult  to  divorce  in  the  measurements. 

e.  Once  the  four  component  losses  mentioned  above  have  a 
firm  basis  in  measurements,  it  will  be  appropriate  to  consider 
Reynolds  number  and  Mach  number  effects  as  they  affect  the  losses  in 
turbine  blade  rows. 

3.   The  computer  program  for  the  turbine  performance  analysis 
yields  results  as  accurate  as  the  loss  data  available.   The  program 
only  requires  inputs  of  loss  coefficients  in  the  specified  form  and 
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a  design  based  on  the  flow  function  CD   for  area  determinations. 

The  analysis  of  the  off-design  performance  of  this  particular 
turbine  indicates  that  large  gains  in  power  output  can  be  obtained 
at  pressure  ratios  two  or  three  times  the  design  value  if  a  157„  drop 
in  efficiency  can  be  accepted.   Attempting  to  obtain  the  last 
possible  increment  of  power,  however,  requires  extreme  values  of 
pressure  ratio  and  another  large  loss  in  efficiency. 

4.   Further  investigation  into  the  after-expansion  behind  a 
blade  row  is  needed.   A  cursory  examination  of  the  magnitude  of  the 
axial  component  of  the  Mach  number  of  the  flow  leaving  an  after- 
expansion  showed  it  ot  be  less  than  one  when  the  loss  coefficient 
decreased  to  zero.   The  exact  reason  for  this  loss  coefficient  be- 
havior is  not  known. 

B.   Recommendations 

1.   That  an  organized  effort  be  made  to  measure  and  define 
component  losses  in  blade  rows  in  a  form  useful  for  design  purposes. 

a.  It  would  appear  that  trailing  edge  thickness  effects 
and  profile  losses  could  be  measured  in  cascade  test  rigs  with  a 
limited  number  of  profiles,  solidities  and  blade  heights.   Some 
correlation  of  these  measurements  could  be  made  with  Ainley  and 
Markov. 

b.  The  tip  clearance  loss  should  be  measured  in  a  rotating 
turbine  test  stand.   The  direct  effects  would  not  be  hard  to  measure 
if  the  effect  on  secondary  loss  was  not  present.   Since  this  effect 
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is  present,  it  will  require  study  to  obtain  quantitative  correlation 
and  isolate  the  tip  clearance  loss  from  the  secondary  flow  effects. 
c.   The  secondary  loss  will  probably  have  to  be  deduced  by 
subtracting  the  outer  component  losses  from  the  overall  efficiency 
as  measured  for  the  complete  turbine  blade  row.   No  method  is  pres- 
ently known  for  measuring  the  secondary  loss  directly. 

2.   Tests  on  the  after-expansion  behind  a  blade  row  should  be 
made.   Physical  measurements  may  shed  some  light  on  the  reason  for 
the  loss  coefficient  reaching  a  peak  and  then  decreasing  to  negative 
values  as  the  expansion  is  increased.   Further  theoretical  study  is 
also  required,  especially  to  consider  three-dimensional  effects. 
The  testing  program  should  include  this  variable. 
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1.2191     1 

.0200 

.0067 

.6800 

•  6800 

1.2338 

.1810 

.0080 

.6600 

.6600 

1.2485 

.3570 

.0095 

.6400 

•  6400 

1.2632     1 

.5500 

•  0112 

.6200 

.6199 

1.2780     1 

.7620 

.0130, 
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•  6000 

1.2928 

.9930 
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•  5800 

.5799 

1.3076     2 
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1.3225    2 
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.0194 
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1.3525    3 

•  1560 
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TIME,   1 
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ANO   36  SECONC 

S 

FIG.  23 

AFTER-EXPANS/ON  AT  V-/37  AND  cx=75K 
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23.0000 
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76 

INLET  ANGLE 
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2C.U000 
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•  a 
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SPACING 
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.9764 

.9507 
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R 
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.7291 

XI 
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X2 
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Yl 

a 

.60  33 
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s 

.6917 
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X5 
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Y3 
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X6 
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Y5 
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Y6 
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.6527 
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a 

.661387 
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a 

.9u34 

.9361 

CG(X7) 

a 

.9433 

.9360 

CG(Y7) 

a 

.6996 

•   .6994 

CG(Y8) 

a 

.675U 

.6756 

I(X-X) 

a  i 

,06302621  . 

06249739 

KY-Y) 

s  « 

,03595965  . 

03399491 

X(MAX)  1 

a 

.6378. 

.6305 

Y(MAX)  1 

a 

.6619 

.6619 

X(MAX)2 

a 
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.•6186 

Y(MAX)2 

a 

.6876 

.6876 
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a 

.U66U 
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TIMEt   0  MINUTES  AND   33  SECONDS 

FOR    ROTOR   DESIGN 
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2.2681 
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.6527 
.637951 

.6527 
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.626808 

.616C3C 
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.9850 

.9292 

.9225 

.9160 

.9290 

.9223 
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.6993 
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.6991 

.6757 
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06143323  . 

C6090044 

03219183  . 
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FIC,.  30 
DIMENSIONS  OF  ROWR  BLADES 
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TURBINE  PERFORMANCE 
INPUT  VALUES  ANO  .CONSTANTS  COMPUTEO 


FLOW   RATE   *      100.0   LBM/SEC           PO   *      1000 
REFERRED   FLOW   RATE   ■  3.821           GAMMA  «    1. 

•    PSIA 
37 

TO   =      1460.    * 

R      «   421.50 

• 

CI    «           3.6195 
C2    ■      316.9026 
C3   «         .000010 
CU    *       136.3049 

• 

EXP1    » 
EXP2   * 
EXP3   » 
RTO      = 

.2701 
3.7027 
2.7027 

38.21 

PHYSICAL    PARAMETERS 

RPM   ■       13000. 

* 

REFERRED   RPM 

1   *               3U0.226 

MEAN   OIA    AT    ROTOR    INLET   *    23.00    IN      STATOR    BLADE 

MEAN   DIA    AT   ROTOR    EXIT      ■    23.00    IN      ROTOR    BLADE 

STATOR  DISCHARGE   ANGLE   ■      75.00   DEG 

ROTOR    INLET    ANGLE                 »      69.60   DEG 

ROTOR    DISCHARGE   ANGLE      *    -71.57   DEG 

1   HEIGHT   *    1.165    IN 
HEIGHT      »    1.U65    IN 

LOSS   COEFFICIENT    THR0U5H   ROTOR    (FLOWS)    ■ 

.1572 

STATOR   EXIT    PLANE 

M(VE)           VE/RTO         PE/PO 
1.00         121.515         .5380 

TE/TO 
•  85  30 

i 

- 

STATOR   DISCHARGE    PLANE 

M(V1)            V1/RT0         P1/P0 
1.06         131.612         .0600 

T1/T0 
.82  75 

ALFA1 

74.78 

ZETA1 
•  0883 

M(Wl)            W1/RT0      PR1/P0 
.80           99.071         .6955 

TR1/T0 
.92  53 

BETA1 
69.60 

• 

IN   ROTOR    ENTRY 

MtWI)            W1/RT0      PR1/P0 
.80           99.071         .6955 

TR1/T0 
.92  53 

BETA1 
69.60 

FLOW   FUNCTION   COMPUTATIONS 

PHI    COMPUTED   «    .58953 
PRESSURE    RATIO    IMPOSED    «      1.512 
ROTOR    EXIT    PLANE 

PHI   CRITICAL   »   .60837 
CRITICAL    PRESSURE   RATIO   «      1.827 

M(WE)            WE/RTO      PE2/P0 
.73           90.952         .4600 

TE2/T0 
•  84  29 

ROTOR    DISCHARGE    PLANE 

■ 

MIW2)           W2/RT0        P2/P0 
•67           84.389        .0600 

T2/T0 
.85  43 

BETA2 
-71.57 

ZETA2 
.2744 

MIV2)           V2/RT0        P2/P0 
.42           53.105         .4600 

T2/T0 
.8543 

ALFA2 
-59. 8U 

RESULTANT 

POWER    AND  EFFICIENCY 

VU1  (FT/SEC)      VU2CFT/SEC)         P2(P'SIA)      P2TJPSIA) 
4852.61              -175U.49                   460.                 519. 

T2             T2T 
1247.       1288. 

OVERALL   PRESSURE  RATIO 

«      2 •  17%      HP  » 

48682. 

EFFICIENCY  «   .723 
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Ji-*  P/AGRAM  {DR  AFTER-EXPANSMN  PROCESS 
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FIG.  36 

FOWER.  COFFF/C/ENT  vs.  TVRB/NE  EFF/C/EMCY 
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Appendix  I.   Development  of  Flow  Function  Relations 

In  Fig.  32  is  shown  the  adiabatic  expansion  process  in  a 
turbine  from  state  1  to  state  2. 

The  fluid  proceeds  from  state  1  to  state  2'  and  the  isentropic 
enthalpy  drop  is  represented  by  process  (1-2!  ).   The  fluid  proceeds 

from  state  2'  to  state  2  and  the  isentropic  enthalpy  drop  is  process 

* 
(2 '-2  ).   For  the  complete  process  the  isentropic  enthalpy  drop  is 

represented  by  process  1-2.  .   However,  the  sura  of  processes  (1-2!  ) 

* 
+  (2-2  )  is  greater  than  process  (1-2.  )  due  to  the  so-called 

X  s 

"reheat  factor".  The  frictional  effects  raise  the  temperature,  and 
therefore  the  enthalpy,  of  the  fluid  as  it  proceeds  from  state  1  to 
state  2'.   This  increase  of  enthalpy  is  then  available  to  do  work. 

As  a  result,  it  is  necessary  to  make  a  distinction  between  the 
loss  coefficient  }  for  the  overall  process  of  Fig.  14  and  the  loss 
coefficient  y       for  the  polytropic  process  from  pressure  (p)  to 
(p-dp)  on  Fig.  32. 

The  polytropic  loss  coefficient  is 


f*  - l  -  ?p  ff.u 


Looking  at  the  expansion  process  in  more  detail  in  Fig.  32, 
the  efficiency  can  also  be  seen  to  be 


fr- 


s£71l  d,2) 
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From  this  the  relation  along  the  polytropic  line  is 


^£  =    * 


'-(*?f 


=  Z 


(1,3) 


and 


<£T  -  ^ 


^7.      = 


yt4t 


Z 


or 


*£T  _    ?„    3W   ^° 


/f 


(1,4) 


Since,  for  an  adiabatic  process 


'7U 


~P/V        »  constant,    then 


r^   ^ 


*^ 


(1,5) 


and  for  an  isentropic  process 


-p<*T        .  constant,  then     /     oC    ^ 


(1,6) 


Comparing  these  exponents,  there  is 


•ru 


-7A^)-c-VM 


and 


STL      = 


^ 


1p(r-i)-+t 


(1,7) 
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To  find  the  relation  between  the  overall  loss  coefficient  and  the 
polytropic  loss  coefficient,  eq.  (1,5)  can  be  written  between  two 
points  on  the  polytropic  expansion  line  as 


51  = 

T 


and 


Z~T=    T0 


*~l 


i-tiL 


3    TZl 


Po 


=  7 


5 


v-i 

-r 


(1,8) 


Solving  equation  (1,8)  for  *h   and  using 


r=  /-? 


the  reduced  result  is 


•f- 


P, 


'  -  ft)* 


(1,9) 


To  define  a  flow  function  for  the  mass  flow  through  the  blade 
row,  the  expression  for  mass  flow  rate  must  be  expressed  in  terms  of 
the  pressure  ratio  and  the  polytropic  exponent.   The  mass  flow  is 


sUS"     — 


AM 


^i 


(1,10) 
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and  following  the  procedure  outlined  in  Appendix  IV  for  finding  the 
velocity  corresponding  to  a  pressure  drop  the  result  is 


r-/ 


/ 


p. 


(IV, 4) 


Along  the  poly tropic  line,  for  an  adiabatic  process 


■? 


■  constant 


or 


When  equations  (IV,4)  and  (1,11)  are  substituted  into  (1,10)  the 
result  can  be  arranged  in  non-dimensional  form  to  define  a  flow 
function  (D   to  be 


t~ 


?o    */// 


r* 


a. id 


(1,12) 


This  flow  function  reaches  a  maximum  at  the  choked  condition  for  the 
blade  row.   To  find  the  corresponding  pressure  ratio  for  choked  flow, 
the  derivative  of  the  flow  function  with  respect  to  the  pressure  ratio 
will  be  set  equal  to  zero  and  the  result  solved  for  the  pressure  ratio. 


f- 


1    = 


_  P 


? 
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Then 


^ 


(**) 


sOl 


zr 


-1-/ 


A  + 


•n 


'ZLtl 


-/ 


A. 


which  reduces  to 


^(0Z)    _  Z* 


^£b 


*-/ 


•* 


/-,* 


When  set  equal  to  zero,  the  result  is 


/t 


z_ 


and 


<&L 


2 

/n+J 


(1,13) 


When  the  relation  for  the  critical  pressure  ratio,  eq.  (1,13), 
is  substituted  into  eq.  (1,12)  the  critical  flow  function  is  obtained 
as 


4   = 


sri-h/ 


Zt  (<>*  -/  \ 


(1,14) 
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Appendix  II.   After-expansion  from  a  Blade  Row 

To  find  relations  to  express  the  expansion  of  a  sonic  or  super- 
sonic flow  out  of  a  blade  row,  the  following  model  is  employed.   The 
cascade  is  assumed  to  have  infinitely  many  blades  of  zero  thickness 
arranged  as  in  Fig.  33  below. 


Fig.  33  Relations  for  After-expansion  from  a  Cascade 

A  unit  height  perpendicular  to  the  page  is  assumed.   Also  the  uniform 
adiabatic  flow  is  assumed  to  be  frictionless  at  (\jj   so  that  R  is 
perpendicular  to  (w)  .    The  control  area  under  consideration  is 
bounded  by  the  surfaces  QtU    ,    (e}   and  (d)  ,  i.e. ,  the  "wall",  the 
"exit"  plane  and  the  "discharge"  plane.   V  is  sonic  or  supersonic. 


-109- 


The  continuity  equation  is  used  to  find  a  relation  between  the 
densities , 

*4-    —     A,  <>&2  <* 


and  so 


The  momentum  equation,  Ref.  1,  ignoring  shearing  stresses  is 


(11,1) 


(11,2) 


where 


A  A 

A  A  A  y 

and  so 

A 

Of  particular  interest  are  the  relations  in  the  i  direction 
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and  using 


^5    =  /'A    JL 


/* 


V.4- 


reduces   to 


/*£ 


^       L_ 


vx—  y^  /ten.  ac< 


=  o 


Since 


v 


/#. 


*MX 


1 


£ 
A 


=    QL  _/ 


(11,3) 


Since  the  process  is  assumed  adiabatic,  the  total  enthalpy  remains 
constant  through  the  control  area, 

z  z 

and 


~JL  =    Ce  / 


=     &       T      -r     =     ^        ^ 


f 


and  substituting 


(11,4) 
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Rearranging  and  substituting  in  eq.  (11,1)  and  using 


■ft  r  =  *J/2. 


results    in 


ft,  =  ^L 


Yjl 


/C&&  o( 


£&!, 


(o(~Ao<) 


V, 


1  +   lttAiZ(/-JU 


Z  n 


-        (11,5) 


Then  equating  eqs.  (11,3)  and  (11,5)  and  collecting  terms  results 
in  a  quadratic  equation  in  velocity  ratio, 


V, 


K 


£ 


VM^CmlAoC-.  ***** 


/O*t(o(-&o() 


V1t+/ 


•/- 


&  4 

(>*%' 


0 

(II. 6) 


The  loss  coefficient  associated  with  the  flow  through  the 
control  area  can  be  defined  to  be 


f  = 

'A 


tf  -  vj 


*e£ 


v 


based  on  the  relations  in  Fig.  34.   Using  perfect  gas  relations  and 
Fig.  34  as  a  reference,  the  following  relations  are  obtained  for  the 
loss  coefficient. 
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£  ik  =  //--^  -  £ 


-^C-^C 


£  = 


^ 


d-^m 


"it 


6  = 


-2     -«, 


\A.X)       \4~s) 


vy 
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/       z 


rs 


•<*, 


^-, 


Z, 


Ik 


'  2.     "A 


7* 


5"  = 


awU 


rtj 


/^i^/« 


6e 


-d-r 


_/*lVfc  u 


£ 


&. 
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z  z  * 
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+    Z    *   K 

j  r-t 


/- 


2k 


Vk  _  £ 


?-) 


f  = 


1 

i 

2       / 

_/*** 

/  + 

Z      1 

& 


^    H 


%^' 


/  " 


r*-  \    yy\ 
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-v 


£-/- 


& 

k 


/    7^ 


-e     / 


r-/  /%- 


/  - 


a  1  f 


(11,8) 
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The  Mach  number  of  flow  behind  after-expansion  is  then,  using  Fig.  34 
as  reference  for  relations, 


Further 


*u 

H- 

«* 

3 

1^ 


4  A*. 


^d/l^U 


(11,9) 
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Substitute  into  (11,9)  above 


1*- 


Vl\Ma 


V, 


-ik 


I- 


(11,10) 


This  computation  is  included  in  the  computer  program  but  is  not 
printed  out.   It  is  inserted  for  ready  use  if  a  check  is  needed. 
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Appendix  III 


MAIN  PROGRAM  FOR  TURBINE  PERFORMANCE  ANALYSIS 
FLOW  CHART 
FORTRAN  NAMES 
PROGRAM  LISTING 

SUBROUTINE  AFTER 
FLOW  CHART 
FORTRAN  NAMES 
PROGRAM  LISTING 

SUBROUTINE  TRNGL 
FLOW  CHART 
FORTRAN  NAMES 
PROGRAM  LISTING 

SUBROUTINE  BEFORE 
FLOW  CHART 
PROGRAM  NAMES 
PROGRAM  LISTING 

SUBROUTINE  CURVE 
FLOW  CHART 
FORTRAN  NAMES 
PROGRAM  LISTING 

SUBROUTINE  CPHI 
FLOW  CHART 
FORTRAN  NAMES 
PROGRAM  LISTING 

PROGRAMS  TURPLOT  AND  TURBINE  (HIGH  HEAD  COEFFICIENT  DATA) 
FORTRAN  NAMES 
PROGRAM  TURPLOT  LISTING 
PROGRAM  TURBINE  LISTING 

PROGRAM  FANNO  3  (FLOW  FUNCTION  AND  ZETA  POLYTROPIC) 
FORTRAN  NAMES 
PROGRAM  LISTING 

PROGRAM  FANNO  4  (CHOKED  FLOW  VALUES  OF  FANNO  3) 
FORTRAN  NAMES 
PROGRAM  LISTING 

PROGRAM  FOR  AFTER -EXPANSION  COMPUTATIONS 
FORTRAN  NAMES  SAME  AS  "AFTER"  ABOVE 
PROGRAM  LISTING 
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rev 


|He::.d   in  data 


PP?;;      —      PP?/     _        A    RPM 

iLfcri     —     AXii  <L*   --—  -- 


1 


Compute  constants  and 
exit  plane  conditions  i 


YES 


Is  p-  >  P, 


o 


Subro 

Slice! 


tine  AFTER 
:  condition 


Subroutine  APT^R 

After- expansion  condition 


V71 


V\        >      0      ? 
'  A 


|  Subroutine    CURVE 


Compute   Y.,  ,    lly1 

t    ' 


Subroutine   TRIIGL 

Compute  vj-  ,  :;.n ,  v^  ,  u1 


I    . 


'Is   Mv    <    1.00    ? 


ILL 


i:o 


Subroutine  BEFORE 


— ^is  f3  > 


0    ? 


NO 


Compute  rotor* entry 


conditions 
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Subroutine   CPHI 


"  m  >  a 


'.^.•routine   CPEI 


1:0 


P2/        \P  /    ci 


Compute   conditions  c,t 
rotor  discliG-r^e  uncIioAed 


X 


Set  %e   =    1  .00 


Compute  ro'tor  discharge 
choked  condition 


Cod  mite   HP 


I 


I, 

__s 

-=^\  All    A  p    ? 


>   ITO 


- 


sfi  In 


V 


A  HPK  ?    ) 


::o 


P2    =    P2  -      A      Po 


.  • 

Subroutine  AFTER 
.  "sfter-exoansion  cor.dit: 

.on 

xcs 


IS      f ,         >        0 


HO 


_y 


<?) 
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PROGRAM  PERFORM 


AE2 

A  0 
e2 

AJ 

J 

ALFAB1 

Bl 

ALFA1 

*i 

ALFA2 

«2 

AME 

M 
e 

AMIS 

M. 
is 

AMV1 

\i 

AMV2 

*V2 

AMW1 

*wi 

AMW2 

"W2 

AMWE 

"we 

AMWEX 


BETA1 

P 

BETA2 

P 

BETAB1 

f 

BETAB2 

h 

CP 


FORTRAN   NAMES,    EQUIVALENT   SYMBOLS  AND  DEFINITIONS 

Area  at   exit  plane  of   rotor 

Mechanical   equivalent  of  heat 

Undeflected  flow  angle  of  absolute  velocity 
leaving  stator 

Flow  angle  of  absolute  velocity  leaving 
stator 

Flow  angle  of  absolute  velocity  leaving 
rotor 

Mach  number  of  flow  at  exit  plane  of  stator 

Mach  number  for  isentropic  expansion 

Mach  number  of  absolute  flow  leaving  stator 

Mach  number  of  absolute  flow  leaving  rotor 

Relative  Mach  number  of  flow  entering  rotor 

Relative  Mach  number  of  flow  leaving  rotor 

Relative  Mach  number  of  flow  at  exit  plane 
of  rotor 

Average  Mach  number  of  flow  at  exit  plane 
of  rotor  for  choked  flow 

Flow  angle  of  relative  flow  entering  rotor 

Flow  angle  of  relative  flow  leaving  rotor 

Relative  flow  angle  for  zero  incidence 
entering  rotor 

Undeflected  relative  flow  angle  leaving 
rotor 

Specific  heat  at  constant  pressure 


Bl 


B2 
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CI 

C2 

C3 

C4 

DM1 

Dml 

DM2 

Dm2 

DRPM 

ARPM 

DTE 

(ATe/To) 

DTF 


DTR 


DTE  2 


DT2 


EFF 


Ar/T 


<TR1-T1> 


DTRAT  (AT/T    ) 

o 


(AT    0/T    ) 
e2      o 


(AT9/T    ) 
z     o 


7 


EXPl 

EXP2 

EXP3 

G 

g 

GAM 

* 

HKR 

*ta* 

HKS 


«h! 


(Constant  #1)    =      "/ R/g 

(Constant  #2)    -       V(2  y  gR)/(^-l) 

(Constant  #3)    =  (*2f -l)/(2  If  gR) 

(Constant  #4)    =       IrFgR 

Mean  diameter  of   blade  at   rotor  entrance 

Mean  diameter  of   blade  at   rotor  exit 

Increment   of  rpm 

Temperature   drop   from  inlet   to  exit   plane 
of  stator 

Final   temperature   difference   from   total    to 
static   temperature  after  rotor 

Temperature  drop  between  relative  total 
temperature  ahead  of  rotor  and  static 
temperature  ahead  of  rotor 

Temperature  drop  from  inlet  to  discharge 
plane  of  stator 

Temperature  drop  from  inlet  to  exit  plane 
of  rotor 

Temperature  drop  from  inlet  to  discharge 
plane  of  rotor 

Efficiency  based  on  total  pressure  after 
rotor 

Exponent  #1  -  (V-1)/1T 

Exponent  #2  ■  T/(T~l) 

Exponent  #3  =  l/( T-l) 

Acceleration  of  gravity 

Specific  heat  ratio 

Blade  height  correction  for  rotor  loss 
coefficient 

Blade  height  correction  for  stator  loss 
coefficient 
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HP 

HR2 

^2 

HS1 

hsi 

NMAX 

PE 

«vv 

PE2 

<pe2/Po> 

PHIC 

A 

PHICR 

0c  r 

PPCR 


(P/P) 


cr 


PO 

P 
o 

PPRAT 

PR1/P2 

PRAT 

pl/pe 

PRATO 

VP2 

PRP 

PR1/P1 

PR1 

PRl/Po 

PRIY 

PRl*/Po 

PI 

Pl/Po 

P2 

P2/Po 

P2T 

P2 

K 
RFR 

W  ■/  T  /P 
/   o  o 

RPM 

rpm/  y  T 

Horsepower 

Rotor  blade  height 

Stator  blade  height 

Number  of  data  points  stored  for  a 
particular  loss  coefficient  curve 

Static  pressure  at  exit  plane  of  stator 

Static  pressure  at  exit  plane  of  rotor 

Flow  function  computed  with  input  *2T  and  } 

Critical  flow  function  for  choked  flow  in 
rotor 

Critical  pressure  ratio  for  choked  flow  in 
rotor 

Total  pressure  ahead  of  stator 

Pressure  ratio  through  rotor 

Pressure  ratio  between  stator  exit  plane 
and  stator  discharge  plane. 

Overall  pressure   ratio 

Relative  pressure  ratio 

Relative  total  pressure  ahead  of  rotor 

Relative  total  pressure  behind  shock  in 
rotor  entry 

Static  pressure  between  stator  and  rotor 

Static  pressure  after  rotor 

Static  pressure  after  rotor 

Gas  constant 

Referred  flow  rate 

Referred  rpm 
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RPMA 

RTO 

/T 

'      o 

TE 

T  /T 
e  o 

TEIS 

T   /T 

TE2 

T  0/T 
e2   o 

TI 

TO 

T 
o 

Tl 

T./T 
1  o 

T2 

VTo 

T2T 

T2 

TP2 

P2 

TRAT 

(Tlls/To> 

TRATR 

(T2ls/TRl> 

TR1 

TD1/T 
Rl  o 

TT2 

T 

T1Y 

T  * 

Ul 

VV^o 

U2 

V  l/  To 

UV 

u1/v1 

VE 

v/*T 

VM 

V  /  V  T 
m  '   0 

VU1 

V  ./  V  T 
ul  /   0 

Actual  value  of  rpm 


Temperature  at  exit  plane  of  stator 

Isentropic  temperature  at  exit  plane  of 
stator 

Temperature  at  exit  plane  of  rotor 

True  interval  between  data  points  on  loss 
coefficient  curves 

Total  temperature  ahead  of  stator 

Static  temperature  between  stator  and  rotor 

Static  temperature  after  rotor 

Static  temperature  after  rotor 

Total  pressure  after  rotor 

Isentropic  temperature  ratio  through  stator 

Isentropic  temperature  ratio  through  rotor 

Relative  total  temperature  ahead  of  rotor 

Total    temperature  after  rotor 

Static  temperature  behind  shock  in  rotor 
entry 

Peripheral  speed  of  blade  at  mean  diameter 
at  rotor  entrance 

Peripheral  speed  of  blade  at  mean  diameter 
at  rotor  discharge 

Abscissa  for  rotor  loss  coefficient  data 

Absolute,  velocity  of  flow  at  exit  plane  of 
stator 

Meridional  component  of  absolute  velocity 

Peripheral  component  of  absolute  velocity 
leaving  stator 
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VU2 

Vu27  V  TQ 

VI 

V  i  To 

V2 

v0/  y  t 

2   *   o 

WDOT 

w 

WE 

W  /  V  T 
e   '   o 

WM 

w  /  y  T 

m  '   o 

QUI 

"ui/yx 

WU2 

wu2/  f*7 

Wl 

V  VTo 

W1Y 

W  */  V  T 

1   '   o 

W2 

V/*T 

YY 

ZET 

1 

ZETAA 

i, 

ZETAB 

^B 

ZETAS 

rs 

ZETAA1 

fal 

ZETAA2 

/a2 

ZETAl 

fl 

ZETA2 


f; 


Peripheral  component  of  absolute  velocity 
leaving  rotor 

Absolute  flow  velocity  leaving  stator 

Absolute  flow  velocity  leaving  rotor 

Mass  flow  rate 

Relative  velocity  at  exit  plane  of  rotor 

Meridional  component  of  relative  velocity 

Peripheral  component  of  relative  velocity 
entering  rotor 

Peripheral  component  of  relative  velocity 
leaving  rotor 

Relative  flow  velocity  ahead  of  rotor 

Relative  flow  velocity  behind  entry  shock 
in  rotor 

Relative  flow  velocity  behind  rotor 

Ordinate  of  loss  coefficient  curves 

Loss  coefficient  for  velocity  through  rotor 
including  after-expansion 

After-expansion  loss  coefficient 

Loss  coefficient  due  to  rotor  entry  shock 

Stator  loss  coefficient 

Loss  coefficient  for  flows  through  stator 

Loss  coefficient  for  flows  through  rotor 

Loss  coefficient  for  velocities  through 
stator 

Loss  coefficient  for  velocities  through 
rotor  without  after-expansion 
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PROGRAM  PERFORM 

DIMENSION  YY(4t53),NMAX<4).TI(U) ' 

COMMON  CI ,C2,C3,C4,GAM,EXP1,EXP2,EXP3,P0,T0,R 

READ  281,P0,T0,GAM,R.WD0T   . 

REAO  282, RPMA, DM1 ,DM2,HS1 ,HR2 , ALFAB1 , BETAB1 ,BETAB2,AE2 

READ  283, AME,  PE.,Z  ETAAl  ,ZETAA2  ,  AMIS 

DO  118  K  *    1,  4  r*>. 

READ  284, J, (NMAX( J),TI (J) ) 

NMAX  =  NMAX(J) 

READ  285  ( YY( J , I )  ,  1= 1 , NMAX ) 
118  CONTINUE 

A  J  =  778.  17  ...  ': 

G  =  32.174 

RTO  =  SQRTF(TO) 

RFR  =  WDOT»RTO/PD 

Cl=SQRTF(R/32. 17* ) 

C2=SQRTF(2.*GAM*3  2.174»R/(GAM-1.) ) 

C3=(GAM-1 .  )/(2.»SAM#32.174»R) 

CU  =  SQRTF(GAM»32.1  7U»R) 

EXP1=(GAM-1 ,)/GAM 

EXP2=(GAM)/(GAM-1 . ) 

EXP3=1./(GAM-1. ) 

CP  =  R»GAM/(AJ*(3AM-1.)) 

PE  =  PE/PO 

TEIS  =  PE*»EXP1- 

DTE  =  (l.-TEIS)«( 1.-ZETAA1) 

TE  -  1.-DTE 

VE  =  C2»SQRTF(DTE ) 

CALL  CURVE  ( 1 ,HS1  , HKS, T I , NMAX ,YY) 

CALL  CURVE  ( 2 , HR2 ,HKR , TI , NMAX, YY) 

DO  200  J  =  1,5 

DRPM  =  1000.«FL04TF( J-1 ) 

RPMA  =  13000. -DR'M 

PRINT  250 

PRINT  251 ,WDOT,PD , TO 

PRINT  252, RFR, GAM, R 
■    PRINT  253,C1,EXP1 

PRINT  25U,C2,EXP2 

PRINT  25b,C3,EXP3 

PRINT  256, C4, RTO 

RPM  =  RPMA/RTO 

PRINT  257 

PRINT    258, RPMA, R?M 

PRINT    259,DM1,HS1 

PRINT  260, DM2, MR? 

PRINT  261,ALFAB1 

PRINT  262,BETAB1 

PRINT  263,BETAB2 

PRINT  26U,ZETAA2 

ALFAB1  =  ALFAB1#3. 1415927/180. 

BETAB1  =  BETAB1«3 . 1415927/180. 

BETAB2  =  BETAB2»3. 1415927/180. 

PRINT  265 

PRINT  266,AME,VE,PE,TE 

L  =  2 

DO  198  I  =  1,30 

PI  =  470. 

PI     =    Pl-FLOATFU) »10. 

PI  =  Pl/PO 
111  CONTINUE 

PRAT=P1/PE 

P2  '    PI 

IF     (Pl-PE)     100,131,100 

100  CALL    AFTER    ( AME.P RAT, ALF AB1 , ALFA1 ,ZET AA ) 
IF     (ZETAA)     116,116,103 

GO    TO    103 

101  CONTINUE 
ALFA1=ALFAB1 
ZETAA=0.0 

103  CONTINUE 

CALL    CURVE     ( 3  ,  AMI S,ZETAS, TI , NMAX,  YY  ) 
ZETA1    -    ZETAS*HKS+ZETAA 
TRAT    =    P1»«EXP1 
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DTRAT  =  ( 1.-TRAT) *(1.-ZETA1 ) 

Tl  =  1. -DTRAT 

VI  =  C2»SQRTF(0nAT) 

AMVl  -  V1/(C4«SQ*TF(T1>) 

IF  (L-3)  119,120,120 
120  PRINT  290 

L  =  1 
1  19  CONTINUE 

PRINT  267 

AL  =  ALFAl«180./3. 1415927 

PRINT  291 ,AMVl,Vl , PI , Tl , AL, ZETAl 

CALL  TRNGL  (  1  , V 1, DM  1 , ALFA  1 ,RPM, VM, VU1 , WU1 ,U1 , Wl ,85 TAl ) 

DTR=C3»Wl*Wl 

TR1=THDTR 

PRP  =  (TR1/T1 )»*EXP2 

PR1  =  P1*PRP 

AMWl  =  Wl/(C4»SQRTF(T1)) 

PRINT  269 

BE  =  BETAl*180./3. 1415927 

PRINT  260, AMWl ,W1 ,PR1,TR1 ,BE 

IF( AMWl-1 .0)104,1  04, 105 

104  Wl=Wl*COSF(BETABl-BETAl) 
PRINT  271 
DTR=C3*Wl«Wl 
TR1=THDTR 

PRP  =  (TR1/T1 )»*EXP2 

PR1  =  P1*PRP 

AMWl  =  Wl/(C4«SQ*TF(T1)) 

BETA1  =  BETAB1 

BE  =  BETAl«180./3. 1415927 

PRINT  268, AMWl, Wl ,PR1,TR1,BE 

GO  TO  106 

105  CONTINUE 

CALL  BEFORE  (AMWl  ,BET AB1 , BETA  1 ,P1 , Wl , TR1 , Tl , Wl Y, P\ 1 Y, T 1 Y, ZETAB ) 
PRINT  270 
Wl  =  W1Y 

AMWl  =  W1/(C4»SQ*TF(T1Y)  ) 
PR1  =  PR1Y 
BETA1  =  BETAB1 
BE  =  BETAl»180./3. 1415927 
PRINT  268, AMWl ,W1  ,PR1 ,TR1 ,BE 
IF  (ZETAB)  114,136,106 
114  PRINT  288 
GO  TO  198 

106  CONTINUE 
UV  =  Ul/Vl 

CALL  CURVE  ( 4 , UV, ZETAR, T I , NMAX, YY) 
ZETA2  =  ZETAR«HKR 
PPRAT=PR1/P2 
TRATR  =  1 ./(PPRAT«»EXP1) 
DTE2  =  (TR1-TRAT**TR1 )»( 1.-ZETAA2) 
.  TE2  '    TR1-DTE2 
WE  =  C2*SQRTF(DTE2) 
AMWE  =  WE/(C4*SCRTF(TE2)  ) 
PRINT  272 

CALL  CPHI  (ZETAAZ , PPRAT, PHIC, PHICR, PPCR) 
PRINT  273, PHIC, PH ICR 
PRINT  27i4, PPRAT, PPCR 
IF  (PPRAT-PPCR)  107,115,108 

107  CONTINUE 

DELTA  =  ABSF(PPRAT-PPCR) 
IF  (DELTA-. 0005)1 15,1 15, 122 
122  PE2  =  P2 
PRINT  27b 

PRINT  266, AMWE, W5,PE2,TE2 
PRINT  276 

DT2=(TR1-TRATR*TR 1 )*( 1.-ZETA2) 
T2  =  TR1-DT2 
W2=C2«SQRTF(DT2) 
AMW2  =  W2/(C4»S0RTF(T2)) 
BETA2=BETAB2   • 
BE2=  BETA2»180./3. 1415927 
PRINT  291,AMW2,W2,P2,T2,BE2,ZETA2 
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CALL  TRNGL  (2 , W2,  DM2,  BETA2,RPM, WM, WU2 , VU2, U2, V2 t ALFA2 ) 

PRINT  277 

AMV2  =  V2/(CU*S<mF(T2)) 

AL2=  ALFA2«1S0./3.1U15927 

PRINT  268,AMV2,V2,P2,T2,AL2, 

GO  TO  113 

108  CONTINUE 
PRINT  286 
L  =  6 

1  10  P2  =  P2  +  .00'01 

PPRAT  =  PR1/P2 

CALL  CPHI  (ZETAA2, PPRAT, PHIC, PHICR.PPCR) 

DELTA  =  ABSF(PPRAT-PPCR) 

IF  (DELTA-. 0005)1 21 ,121,  110 
121  PI  =  P2 

GO  TO  111 

115  CONTINUE 
AMWE  •=  1.00 
PE2  =  P2 

PPRAT  =  PR1/PE2 
TRATR  =  1 ./{PPRAT»«EXP1 ) 
DTE2  =  { rRl-TRAT^»TRl)»( 1.-ZETAA2) 
TE2  =  TR1-DTE2 
WE  =  C2»SQRTF(DTE2) 
AMWEX=  WE/(C4*SQRTF(TE2)  ) 
299  FORMAT  (100X.F10.3) 
PRINT  275 

PRINT  266, AMWE, WE ,PE2,TE2 
DT2=(TR1-TRATR»T^1 )*( 1.-ZETA2) 
T2  =  TR1-DT2 
W2=C2*SQRTF(DT2) 
BETA2=BErAR2 

CALL  TRNGL  ( 2 , W2,  DM2, BETA2.RPM, WM, WU2 , VU2 ,U2, V2, AL FA2 ) 
AMW2  =  W2/(CU*SQ*TF(T2)) 
AMV2  =  V2/(CU#SQ*TF(T2>) 
PRINT  276 

BE2=  BETA2»180./3. 1415927 
PRINT  291,AMW2,W2 , P2, T2, BE2, ZETA2 
PRINT  277 

AL2=  ALFA2»180./3. 1415927 
PRINT  268,AMV2,V2,P2,T2,AL2 
L  =  3 
GO  TO  113 

112  P2=P2-.01 

IF     (L-5)     109,123,  123' 
123    PRINT    290 
L    =    1 

109  CONTINUE 
PRINT    287 
PRAT=P2/PE2 

CALL    AFTER    ( AMWE. PRAT , BET AB2, BETA2, ZETAA ) 
IF     (ZETAA)     116,116,117 

116  PRINT    289 
GO    TO    199 

117  CONTINUE 
PPRAT    =    PR1/P2 

TRATR    =    1 ./(PPRAT«*EXP1 ) 

DT2  =  (TRl-TRATR«n 1  )»(  1 .- ( ZETA2  +  ZETAA) ) 

T2  =  TR1-DT2 

W2=C2*SQRTF(DT2) 

CALL  TRNGL  ( 2 , W2,  DM2, BET A2, RPM, WM, WU2 ,VU2,U2, V2, ALFA2 ) 

AMW2  ~    W2/(C4*SQ*TF(T2)) 

AMV2  =  V2/(Ct**SQ*TF(T2)) 

PRINT  276 

BE2=  BETA2»180./3. 1415927 

ZET  =  ZETA2+ZETAA 

PRINT  291 ,AMW2,W2 , P2, T2,  BE2,  ZET 

PRINT  277 

AL2=  ALFA2«180./3. 1415927 

PRINT  268,AMV2,V2,P2,T2, AL2 

113  CONTINUE 
P2T  =  P2*P0 
T2T  =  T2*T0 
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DTF    =    C3*V2»V2 

TT2    =    T2+DTF 

TP2    =    P2T*( (TT2/r2)*»EXP2) 

PRATO    =    1./P2 

TT2    *    TT2*T0 

OH1    =    T0-TT2 

DH2    =    TOM1  .-(  (TP2/P0)**EXP1  )  )  •  ' 

EFF    =    DH1/DH2 

VU1T    =    VUl»RTO 

VU2T    =    VU2»RTO 

HP    =    RTO»WDOT*l .+14»(VU1T»U1-VU2T«U2)/(G*AJ) 

PRINT    270 

PRINT    279,VU1T, VJ2T , P2T , TP2, T2T,  TT2 

PRINT    280, PRATO, HP, EFF 

L    =    L  +  l 

IF  (AMWE-1.00)  198,112,112 
2500FORMAT  ( 1  Hi ,2 IX 19HTURBI NE  PERFDRMANCE//1 3X35H INPUT  VALUES  AND  CONS 

1TANTS  COMPUTED    /) 
2510F0RMAT( 13H  FLOW  RATE  =  F6.1,  17H  LBM/SEC     PO  =  F6.0,  1 3H  PSIA. 

1T0  =  F6.0,  2H  R   ) 
2520FORMAT  (22H  REFERRED  FLOW  RATE  ■  F5. 3 ,4X8HGAMMA  =  F4.2.7X5HR  - 

1  F6.2  /) 

253  FORMAT  {7X5HC1  =  F9.4, 2 1 X7HEXP1  =  F6.4) 

254  FORMAT  (7X5HC2  =  F9. 4 , 21 X7HEXP2  =  F6.4) 

255  FORMAT  ( 7X5HC3  =  F9. 6 , 2 1 X7HEXP3  =  F6.4) 

256  FORMAT  (7X5HCU  =  F9.4, 2 1 X7HRT0   =  F6.2) 

257  F0RMATI/20X20H  PHYSICAL  PARAMETERS  /) 

258  F0RMAT(7H.RPM  =  = 7. 0 , 23X 1 5HREFERRED  RPM  =  F12.3/) 
2590F0RMAT(27H  MEAN  DIA  AT  ROTOR  INLET  =    F5.2,27H  IN   STATOR  BLADE  HE  I 

1GHT  =  F5.3,3H  IN   ) 
2600F0RMAT(27H  MEAN  OIA  AT  ROTOR  EXIT   =  F5.2,27H  IN   ROTOR  BLADE  HEIG 
1HT   =  F5.3.3H  IN  ) 

261  FORMAT{ 12X26H  ST* TOR  DISCHARGE  ANGLE  =  F6.2,4H  DE3  ) 

262  FORMAT( 12X26H  ROT  OR  INLET  ANGLE       =  F6.2.4H  DES  )  - 

263  FORMAT{ 12X26H  ROTOR  DISCHARGE  ANGLE   *  F6.2.4H  DEG/) 

264  FORMAT  (  7X41HL0SS  COEFFICIENT  THROUGH  ROTOR  (FLOdS)  =  F7.4//) 

265  F0RMAT(18H  STATOR  EXIT  PLANE//32H  M(VE)     VE/RTO    PE/PO    TE/TO) 

266  FORMAT  ( F6. 2, F 1 0. 3 , 2F8. 4/ ) 

2670F0RMAT  (23H  STATOR  DISCHARGE  PLANE//49H  M(V1)     /1/RTO    Pl/PO 
1T1/T0    ALFAl    ZETA1   ) 

268  FORMAT  ( F6  .  2, F 1 0. 3, 2F8.4 , F8 . 2/ ) 

269  FORMAT  {  41H  M(Wl)     Wl/RTO   PR1/P0   TR1/T0    BETA1   .) 
2700F0RMAT(29H  IN  ROTOR  BEHIND  ENTRY  SH0CK//41H  M(Wl)     Wl/RTO   PRl/P 

10   TR1/T0    BETA1    ) 
2710F0RMAT( 15H  IN  ROTOR  ENTRY//41H  M(Wl)     Wl/RTO   PR1/P0   TR1/T0    B 
1ETA1   ) 

272  F0RMAT(27H  FLOW  FUNCTION  COMPUTATIONS) 

273  FORMAT  (10X16H  PHI  COMPUTED  =  F6. 5 , 1 3X15HPHI  CRITICAL  =  F6.5) 
2740F0RMATI26H  PRESSJRE  RATIO  IMPOSED  =  F6.3,28H   CRITICAL  PRESSURE  RA 

1TI0  =  F6.3) 
275  FORMAT  ( 1 7H  ROTOR  EXIT  PLANE//32H  M(WE)     WE/RTO   PE2/P0   TE2/T0) 
2760FOPMAT  (22H  ROTOR  DISCHARGE  PLANE//49H  M(W2)     W2/RT0    P2/P0    T 

12/TO    BETA2    ZETA2   ) 
277  F0RMATI41H  M(V2)     V2/RT0    P2/P0    T2/T0    ALFA2   ) 
2780F0RMAT  (  /15X31H  RESULTANT  POWER  AND  EFFICIENCY  // 

3X1  1HVUKFT/SE:  ),2XHHVU2(FT/SEC),3X8HP2(PSIA),2X9HP2T(PSIA),3X 

2  2HT2,5X3HT2T) 

279  FORMAT  {  F  1  2  .2 , F 1 3 . 2 ,F 1 1 .0 , 2F 10. 0, F7.0/ ) 

280  FORMAT  ( 5X25HCVER ALL  PRESSURE  RATIO  '    F6.3.2X5HHP  =  F6.0,15H   EFFI 
1CIENCY  =  F4.3///) 

281  FORMAT  (5F10.0) 

282  FORMAT  (9F8.0) 

283  FORMAT  (5F10.0) 

284  FORMAT{ 13,16, F 1 0.  0)   ' 

285  FORMAT  (6F9.0) 

266  F0RMAT(/48H  ROTOR  CANNOT  SATISFY  FLOW  FUNCTION,  INCREASE  Pi  /) 

287  FORMAT  {21H  ROTOR  AFTEREXPANS ION) 

288  FORMAT  ( 1 9H  ZETAB  IS  NEGATIVE.   ) 

289  FORMAT  ( 1 9H  ZETAA  IS  NEGATIVE   ) 

290  FORMAT  ( lHl ) 

291  FORMAT  ( F6. 2, F 10. 3, 2F8. 4, F8. 2 ,F8.4/ ) 

198  CONTINUE 

199  CONTINUE 
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200 


ALFAB1  = 
BETAB1  * 
BETAB2  = 
CONTINUE 
END 


ALFAB1*! 80./3.1U15927 
BETAB1»180./3.1U15927 
BETAB2*! 80./3.1U15927 
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I.  PUT 


■c-vr/^, 


14         c 


r 


• 

s 

=.(  -.1)  "■' 

...  ,  -.:     0 

' 

' 

1 

I 
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SUBROUTINE  AFTER 


ZA 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 
Undeflected  flow  angjgfi  (blade  angle) 
Cosine  ratio  of  blade  angle  to  flow  angle 
C^B  -  0(       (degrees) 

°^B  "  °^1   (radius) 
Imposed  pressure  ratio 

Computed  pressure  ratio 

Velocity  at  discharge  plane 

Velocity  ratio  of  flow  at  discharge  plane 
to  exit  plane 

Loss  coefficient  due  to  after-expansion 


ALFA 

*B 

COSRAT 

DALFA 

Ao< 

DALFAR 

Ao< 

PRAT 

(Pd/Pe) 

PRAT  I 

<Pd<Vi 

VD 

vd 

VELRAT 

Vd/Ve 

i. 
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uoo 

401 


402 


410 
409 
411 


403 
404 
405 


406 


412 
413 


407 
408 

414 

415 


SU 
CO 
AL 
£E 
DA 
M 
M 

DA 
DA 
CO 
AA 
BB 
CC 
DD 
IF 
VE 
OPR 
1  / 
IF 
PR 
FO 
CO 
ZA 
VD 
DI 
IF 
IF 
M  = 
DA 
DA 
CO 
AA 
BB 
CC 
DD 
IF 
VE 
OPR 
1  / 
IF 
PR 
CO 
ZA 
VD 
DI 
IF 
IF 
CO 
IF 
AL 
RE 
AL 
ENI 


HROUTIN 
MMON  CI 
FA  =  AB 
=  (2./( 
LFA  =  0 
■  0 
=  M  +  l 
LFA  =  ( 
LFAR  = 
SRAT  = 
=  GAM* 
=  -GAM 
=  COSR 
=  BB*B 
(DD)  4 
LRAT  = 
ATI  =  ( 
VELRAT 

(PRAT1 

INT  409 

RMAT  (7 

NTINUE 

=  l.-{ 

-  VE*V 

FF  =  AB 

(DIFF- 

(PRAT1 

M+l 

LFA  =  { 
LFAR  = 
SRAT  = 
=  GAM* 
=  -GAM 
=  COSR 
=  BB*B 
(DD)  40 
LRAT  = 
ATI  =  ( 
VELRAT 

(PRAT1 

INT  409 

NTINUE 

=  l.-( 

=  VE*V 

FF  =  AB 

(DIFF- 

(PRAT-P 

NTINUE 

(ALFAS 
FA1  =  A 
TURN 
FA1  =  A 
D 


E  AFTER  (AME, PRAT, ALFAS, ALFA1,ZA) 

,C2,C3,;4,GAM,EXP1,EXP2,EXP3,P0,T0,R 

SFtALFAS ) 

(GAM-1.) *AME*AME) ) 

.0 

-.1)**M-1)  +  DALFA 

DALFA*3. 1415927/180. 

C0SF(ALFA)/C0SF( ALFA-DALFAR) 

AME* AME»CCSF(DALFAR) -COSR ATM  GAM- l.)*AME*AME/2. 

*AME*AME-1. 

AT*( l.+( (GAM-1. )*AME*AME/2.) ) 

B-4.*AA»CC 

01,402    , 402 

(-BB+SQ*TF(DD) ) /(2.*AA) 

C0SRAT*(  l.  +  ( (GAM-1.)*AME*AME/2.)*(1.-VELRAT*VELRAT)  ) ) 


)     410,41 

,N,AME,P 
2X,I3,4F 


1,411 

RAT,ALFA,DALFAR 

10.5) 


VELRAT*^ELRAT)/(1.+EE-EE*(PRAT1**EXP1) ) 

ELRAT 
SF(PRAT1 
.00001 ) 
-PRAT)  '4 


-PRAT) 

408,408.403 

04,408,401 


-.  1  )**M-1  )+DALFA 

DALFA*3. 1415927/180. 

COSF ( ALF A )/C0SF( ALFA-DALFAR) 

AME*AME»C0SF(DALFAR)-C0SRAT*(GAM-l.)*AME»AME/2. 

*AME*AME-1. 

AT»( l.-K  (GAM-1. )*AME*AME/2.) ) 

B-4.*AA»CC 

0,406,436 

(-BB+SQ*TF(DD) )/(2.*AA) 

COSR AT* (  l.+  ( ( GAM-1.)* AME* AME/2.  ) * ( 1 .-VELR AT* VELRAT )  ) ) 

)    412,413,413 
,N,AME,?RAT,ALFA,DALFAR 


VELR  AT  *«/ ELRAT) /(  1  .  +  EE-EE*  (  PR  AT  1  **EXP  1 )  ) 

ELRAT 

SFIPRAT1 

.00001 ) 

RATI)     43 


-PRAT) 

408,408,407 

0,408,405 


)     414,41 
LFAS+DAL 


5,415 
FAR 


LFAS    -    OALFAR 
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--    T) --r-! 

1  ,    T,    , 

1    OUTPUT 

t 

1 

'               ■•''.■■ 

31=1     ib solute  values  input 

.'1-id.  relative 
17  =  2       1    fcive   values  input 

. ind  absolute 

' 

r    • 

• 

Gom  rut.e            / 

U.    7        V          P       V 
■'      r'      Ui'*  r .  '      1 

* 

• 

' 

' 

1 

/JTURI? 
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B 

N 

U 

VM 

V 
m 

VU1 

Vul 

WU1 

Wul 

HI 

wl 

SUBROUTINE   TRNGL 
FORTRAN  NAMES,   EQUIVALENT   SYMBOLS  AND  DEFINITIONS 

Input  angle  of   flow 

Output   flow  angle    (If  A  ■  o(    ,    B  =  tf      and   if  A  =/ 
B   -*\)  r  ' 

Value  selects  plus  or  minus  SIGN 

Peripheral  velocity  of  blade  at  mean  diameter 

Meridional  component  of  absolute  velocity 

Peripheral  component  of  absolute  velocity 

Peripheral  component  of  relative  velocity 

Relative  velocity 
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SUBROUTINE  TRNGL  ( N, V 1 ,DM, A, RPM. VM, VUl , WUl ,U, Wl , B) 

COMMON  CI ,C2,C3,:u,GAM,EXP1f EXP2,EXP3,PO,TO,R 

SIGN  =  (-1  )»*N 

U  =  3.1415927*RPM«DM/720. 

VM  -  VUCOSF(A) 

IF  (A)  701,702,732 

701  VUl  =  -Vl«SINF(A) 
WUl  =  VUl  +SIGN*U 
E  =  ATANF(WU1/VM) 
Wl  =  WU1/SINF(B) 

702  CONTINUE 

VUl  =  V1«SINF(A) 
WUl  =  VU1+SIGN»U 
B  =  ATANF(WU1/VM) 
Wl  =  WU1/SINF(B) 
END 
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:::::;ut 
fcjach,  oC 

,<*  ,  v  p,. 

OUTPUT... 

m          nv  *          "D     *'" 

i 

*iH 

•  °3 
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SUBROUTINE  BEFORE 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


A 

ALFA 

ALFAI 

AMD 

B 

CC 

D 

DALFA 
DTR 

FF 

GG 


*i 


M 


&o( 


(Trl-Tl) 


pp 

P1/P1* 

PRP 

PR1/P1 

PR1Y 

PRl*/Po 

P1Y 

Pl*/Po 

W 

w1*/w1 

Interim  computation  for  velocity  ratio 

Blade  angle  entering  rotor 

Flow  angle  ahead  of  rotor 

Mach  number  of  flow  after  shock 

Interim  computation  for  velocity  ratio 

Cosine  ratio  of  angles  involved 

2 
A  -  B,  used  to  determine  whether  shock  should 

be  computed,  or  the  component  of  the  velocity 

in  direction  of  blade  be  used. 


Ac*   =  oC1  -  o(} 


ZETAB 


*. 


Temperature  difference  between  relative  total 
temperature  and  static  temperature 

Interim  computation  to  find  Mach  number  after 
shock 

Interim  computation  to  find  Mach  number  after 
shock 

Static  pressure  ratio  through  shock 

Pressure  ratio  through  shock 

Relative  total  pressure  behind  shock 

Static  pressure  behind  shock 

Velocity  ratio  of  flow  ahead  of  and  behind 
shock 

Loss  coefficient  of  flow  due  to  entry  shock 
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SUBROUTINE  BEFORE  ( AM  1 , ALFA , ALFA1 , PI , Wl , TR1 , T 1 , W 1Y , PR1 Y, Tl Y, ZETAB) 

COMMON  CI  ,C2,C3,C4,GAM,EXP1,  EXP2, EXP3 ,P0 ,TO,R 

OALFA    =    ALFA1-ALFA 

C:    =    COSF (ALFA! )/COSF(ALFA) 

A    «     (  l./(GAM+l.  ))*[  (  1./(CC»AM1»AM1) ) +GAM*COSF { DALF A ) ) 

B    =    -U.-(2'./(GAM  +  1.  )  )*(  (  l./(  AM1»AM1)  )+GAM)) 

D    =    A«A-R 

IF    (D)    602. 60 1.. 631 

W1Y    =    Wl*COSF( ALFA1-ALFA) 

=  0.0  • 

C3»W1Y*W1Y 

TR1-DTR 

<TR1/T1Y)«»EXP2 
P1Y*PRP 


ZETAB 

DTR  = 

T1Y  = 

P*P  = 

PR1Y  = 

RETURN 

CONTINUE 

VV  =  A-SORTF( A*A-B) 

PP=CC»( ( 1 .+( GAM-1 •) «AM1*AM1/2. )•(!•/ VV)-( (GAM-1 


FF  =  (GAM-1 .)«AM1 »AMl/2. 

GG  =  l.-VV»VV 

AMD  =  VV*AMl*SQRrF(1./(1 .+FF«GG) ) 

ZETAB  =  1 .-VV»VV- (2./( (GAM-1 . )*AMl*AMl ) )»( ( (PP)»»EXP1 )-l 

WlY  =  W1*VV 

P1Y  =  P1«PP 

T1Y  =    TRl-C3*WlY*WlY 

PRP  =  (TR1/T1Y)»»EXP2 

PR1Y  =  P1Y«PRP 

END 


J»AMl*AMl»VV/2.) ) 
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-     ^c\v:t:'}L  o':^rE 


j  OUTPUT 


JT 
1:,    JUSsc?issa.{  Ordinate 


Cell   In  array    (at) 


Write  equation  of  parabola 
through  three  nearest  points 


Commute  ordinate 


RSTUEM 


\! 
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SUBROUTINE  CURVE 
FORTRAN  NAMES,.  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


AA 

BB 

CC 

1 

J 

K 

M 

M 

NMAX 

ORD 

Tl 
12 

X 

XA 


TI/2 


The  terms  of  equation  ORD  ■  a  +  bx  +  ex 


Subscripts  tor  ordinate  data  points 


Counter  to  locate  three  nearest  data  points 
to  input  abscissa 

Array  selector 

Number  of  data  points      iected  array 

Ordinate  computed,  i.e.,  loss  coefficient 
value  desired 

Data  abscissa  interval 

One-half  of  data  interval.   Used  to  select  the 
three  nearest  points 

Abscissa  in  terms  of  three  nearest  points 

Input  abscissa 
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SUBROUTINE  CURVEt  N,XAl ,ORD, T I , NMAX, YY ) 
DIMENSION  YY(U,50  )  ,NMAX(4).TlU)  ,Y(50) 
COMMON  C1,C2,C3,:4,GAM, EXP1 , EXP2, EXP3, PO ,TO, R 
XA  =  XA1 
NMAX=NMAX(N)   ■ 
TI=TI(N) 

00  306  I  =  1,NMAX 
306  Y( I )=  YY(N, I) 
M  =  XA/TI 
X  =  T  I*FL0ATF(M) 
M  =  M  +  l 
IF  (M-l)  304,300,301 

304  PRINT  212 
RETURN 

300  I  =  M 

J  =  M+l 
K  =  M+2 
GO  TO  305 

301  NM  =  NMAX-1 

IF  (  M-NM)  302,50  3,303 
303  I  =  M-l 
J  «  M 
K  =  M  +  l 
GD  TO  305 

302  OIFR  =  XA-X 
T2  =  TI/2. 

IF  (0IFR-T2)  303,  303,300 

305  CONTINUE 
AA  =  Y( I) 

CC  =  0.5»(Y(K)+Y(  I )-2.»Y( J)  ) 
BB  =  Y( J)-Y(I  )-C; 
XA  =  (XA/TI)-FL0ATF(I-1) 
ORD  =  AA+BB«XA+C:*XA*XA 
212  FORMAT  (36H  ABSCISSA  WRONG  IN  SUBROUTINE  CURVE 
END 
END 
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£.p/p 


•       OjJTPUT 


Os. 


ft-< 


■D    +    .001 


Compute  ^ 


NO 


is  ?>  ^  i  y 


YES 


Compute  >^ 


HO 


Is    *><,      ? 


J .  i~j 


Compute    ¥  ,    for 

■  (*/P>.or 

1 

' 

RETURN 
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SUBROUTINE  CPHI 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


A 

CPHI 

DAB 

DAN 

DIFF 

EN 

n 

EXP1 

EXP4 

EXP6 

EXP7 

GAM 

T 

PHI 

4> 

PHI  2 

<? 

PHICR 

4 

r  cr 

PPCR 

(P/P) 

RC 

RA 

RB 

ZETA 

f 

ZETAP 

& 

Interim  critical  pressure  ratio  while  computing 
Compute  (p ,    subroutine  name 


cr 


Name   for   the    term   (p/P) 
Name   for   the   term   (p/P) 


exp7 
exp6 


Difference  between  given  and  computed  loss 
coefficient 

Polytropic  exponent 

Exponent  #1  -  (  V   -  1)/  W 

Exponent  #4  -  (1  -  J  )  (  V  -  1)/^ 

Exponent  #6  =  2/n 

Exponent  #7  ■  (n+l)/n 

Specific  heat  ratio 

Flow  function  computed 

Flow  function,  squared 

Choked  flow  value  of  flow  function 

Choked  flow  value  of  pressure  ratio 

Computed  loss  coefficient 

exD^ 
Name  for  the  term  (p/P) 

Name  for  the  term  (p/P) 

Loss  coefficient 

Polytropic  loss  coefficient 
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SUBROUTINE    CPHI    ( ZETA, A, PHI , PHICR, PPCR) 
COMMON   C1 ,C2,C3,:U,GAM,EXP1,EXP2,EXP3,P0,T0,R 
ZETAP    «    -.001 
502    ZETAP    *    .OOUZETAP 

EXPU    -    (1. -ZETAP)  MGAM-1  .  )/GAM 
RA    -     1  ./(  A**FXPM 

Ar   =    i./(a«»;;xpi)   . 

RO    (RA-RB)/(  l.-*B) 

01 FF    =    ZETA-RC 

IF    (OIFF)    501,503,502 
501    ZETAP    *    ZETAP-. 03001 

EXPU  =  (1 .-ZETAP) »<GAM-1.)/GAM 

RA  *  1./(A*»EXP4) 

RB  =  1./<A**EXP1) 

RC=  (RA-RB)/( l.-*B) 

OIFF  *  ZETA-RC 

IF    (OIFF)    501,503 ,500 
500    EN    *    GAM/( 1 .+ZET4PM GAM-1 . ) ) 

EXP6  *  2. /EN 

EXP7  *  (EN«-1.  )/EM 

DAN  *  l./(  A«»EXP6.) 

DAB  »  l./(A*»EXP7 ) 

PHI2  *  2.»(GAM/(SAM-1.))«(0AN-0AB) 

PHI  *  S0RTF(PHI2) 

PHICR  *  (  (2./(EN«-l.))»*(  l./(EN-1.)))»SQRTF(2.»GAM»(EN-l.)/ 
1   < (GAM-1.)»(EN+1.))) 

PPCR  *  {(EN*1.)/2.)»»(EN/CEN-1.)) 

END 
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PROGRAMS  TURPLOT  AND  TURBINE 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


AKE 

AKIS 

ALFA1 

ALFA2 


k 


is 


<*. 


BETA1 

/i 

BETA2 

/2 

DBETA 

¥ 

DETA 

*1 

EFF 

ETA 

7 

ETAA 

h 

ETAC 

% 

EXX 

X 

NALFA1 

«i 

NPHIE 

A, 

PHI 

f 

PHIR 

A 

Leaving  loss  coefficient 

Head  coefficient,  eq.  1 

Angle  of  flow  exiting  stator  (radians), 
absolute 

Angle  of  flow  leaving  rotor  (radians) , 
absolute 

Flow  angle  entering  rotor  (degrees), 
relative 

Flow  angle  leaving  rotor  (degrees)  , 
relative 


PSI 


¥  "A -fit 

n  - 1  -  ?a 

Reheat  factor,  defined  in  Ref.  1 

Efficiency 

Efficiency  assumed 

Efficiency  computed  as  an  interim  step  for 
comparison 

Expression  defined  by  eq.  5a 

Angle  of  flow  exiting  rotor  (degrees), 
absolute 

Carryover  coefficient  of  kinetic  energy  into 
next  stage 

Velocity  coefficient  in  stator 

Carryover  coefficient  of  kinetic  energy  from 
stator  into  rotor 

Velocity  coefficient  in  rotor 
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PSIA 


RADRAT 

R2/R1 

RBETA1 

A 

RBETA2 

A 

RELRAT 

w  /w 

V"l 

RSTAR 

r* 

VELM1 

Vml/Ul 

VELRAT 

<Vm2/V, 

WSQR1 

W  2 
"l 

WSQR2 

W  2 
W2 

XSQR 

X2 

Velocity  coefficient  in  rotor  computed  for 
comparison 

Radius  ratio  through  rotor 

Flow  angle  entering  rotor  (radians),  relative 

Flow  angle  leaving  rotor  (radians),  relative 

Ratio  of  relative  velocities 

Degree  of  reaction,  that  fraction  of  the 
expansion  through  the  stage  which  occurs  in 
the  rotor 

Meridional  component  of  absolute  velocity 
leaving  stator/Peripheral  velocity  of  rotor 
at  mean  diameter 

Velocity  ratio  of  meridional  components 

The  relative  velocity  entering  rotor,  squared 

The  relative  velocity  leaving  rotor,  squared 

Expression  defined  by  eq.  5a 
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PROGRAM  TURPLOT 

DO  50  NALFAl  =  60,75,5  ' 

=  1 
DO  50  I  «  1,12  ' 

RSTAR  —  . 05*(FLOATF(I)  )•. 05 
IF(N-2)21,22,22  ■  • 

21  PRINT  UO 
N  *  N+1 
GO  TO  23 

22  M  -  1 
PRINT  U7 

23  CONTINUE 
NPHIE  *  0 
RADRAT  *  1.00 
VELRAT  *  1.00 
PHI  '    0.95 
PHI R  =  0.95 
EFF  =  O.C 
ETAC  ■  .UO 
PSI   *  .99 

ALFA1  =  FL0ATF(NMFA1  )«3. 1 Ul 5927/ 1 80. 

AKIS    *    2.0 

J    =    1 
13    AKIS    *    AKIS    ■*■    2.0 

J    =    J    ♦    1  . 

170RBETA1    *    ATANF(TANFULFAT)    -    (( 1 .0) /( PHI »COSF { ALFA 1)*SQRTF 

1  (1.    -    RSTAR)»S3RTF(AKIS)  )  )) 
BETA1    *    RBETAU130./3.1U15927 

12    CONTINUE 

OXSGR  *  (RSTAR*( l.+EFF))  ♦  (( PHI*PHI )•( 1 .-RSTAR)»( PHIR-t VELRAT* 
VELRAT«C0SF(ALFA1)»C0SF(ALFA1)/(PSI»PSI))))  -  (2.»PHIR«PHI 

2  »S  INF ( ALFA 1)»SQ*TF( 1 .-RSTAR) /SQRTF ( AKIS ))♦(( PHI R*RADRAT» 

3  RADRAT-1. )/AKIS ) 

19  EXX    =    SQRTF    (  XSQ* ) 
RPETA2=ATANF(-PSI»EXX/(PHI«VELRAT#C0SF(ALFA1)«SQRTF(1.-RSTAR)) ) 
BETA2    =    RBETA2»130./3.1U15927 

20  DBETA    =    BETA1-BETA2 

PSIA  ■■  0.99  -  .000223*DBETA  -  4.97/ ( 1 80. -DBETA) 

DPSI  *  ABSFtPSI  -PSIA) 

IF  (.0000001-DPSI )  10,11,11 

10  PSI  -  PSIA 
GO  TO  12 

11  PSI  =  PSIA 

OETA  =  (2./SQRTF(AKIS) )*(  ( PHI «SINF ( ALFA1 ) »SQRTF( 1. -RSTAR) )- 

1  (RADRAT»RADRAT/SORTF( AKIS) ) ♦ ( PS I*EXX»RADRAT ) ) 

OAKE  »AKIS»( (PHI»»HI»( 1 .-RSTAR ) *CDSF ( ALFA  1) *COSF ( ALFA1 ) ♦VELRAT* 
VELRAT)-M(PSI»PSI)*XSQR)-(2.*PSI»EXX*RADRAT/SQ*TF(AKIS)  ) 

2  *•( (RADRAT*RAD*AT)/AKIS) ) 

VELM1  =  PHI»SORT= ( 1 . -RSTAR ) »SQRTF (AKI S) *COSF( ALFA1 ) 
0ALFA2  -  ATANF( ( R AORAT/SQRTF( AK IS ) )-( PSI »EXX) )/ ( PH I»VELRAT* 
1    COSFIALFAl )»S3RTF( 1. -RSTAR))) 

ALFA2  *  ALFA2*13D./3.1U15927 
0WSQR1  =  1.>(AKIS»PHI*PHI*(1.-RSTAR))-2.*SQRTF(AKIS )»PHI* 
1    SQRTF( l.-RSTA* )*SINF( ALFAl ) 
0WSQR2  =  (PSI*PSI»AKIS)»( ( (1 ,+EFF )«RSTAR) *(PHIR»PHI «PHI» 

1  (1 .-RSTAR) )-(2.»PHIR*PHI»SQRTF(l.-RSTAR)«SINF(ALFAl ) )/ 

2  (SQRTF  (AKIS))*-  ( PHIR+R ADRAT«RADRAT-1 . )/AKIS ) 
RELRAT  »  SQRTF(W5QR2)/SQRTF(WSQR1 ) 

IF  (J-6)  14,14,15 
15  PRINT  U7 

J  =  2 
140 PR  I  NT  41 , NALFA1,MPHIE,RADRAT, VELRAT, PHI, PHI R, PSI, *  STAR, AKIS, 
1 ET A, AKE,VELM1, RELRAT, BET Al, BET A2,ALFA2 
IF  (ETA-ETAC)  50,50,13 
400F0RMAT(/118H1  Al    PE    R2/R1  ■   VM2/VM1    PHI    PHIR     PSI     R* 

1     KIS       ETA      KE    VM1/U1   W2/W1    BETAl    BETA2   ALFA2   /) 
410F0RMAT  (  I3,I6,F7.1,F9.1,F8.2,F7.2>2F7.3,F8.3,F9.5,2F7.3, 

1  F8.3,2F8.2,F7.2  ) 
U7  FORMAT  (/) 
50  CONTINUE 
END 
END 
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PROGRAM  TURBINE 

00  50  N«l,4 

NALFA1  «  80  -  (5»N) 

PRINT  UO 

DO  50  1*1,. 12 

RSTAR*-. 05+(FLOATF( I) )».05 

IF  (1-7)  13,20,13 
20  PRINT  UO 
18  CONTINUE 

J  -  1 

K  =  1 

NPHIE  *  0 

ALFA1  =  FL0ATF(NALFA1)»3. 1415927/180. 

RADRAT  =  1.00 

VELRAT  =  1.00 

PHI  =  0.95 

PHIR  *  0.95 

AKIS«8.0 

EFF  *  0.0 

PSI   *  .99 

ETAA  =  0.0 
170RBETA1  *  ATANF(TANF(ALFA1 )  -  ((  1.0)  /  (PHI«COSF( ALFAl ) »SQRTF 

1  (1.    -    RSTAR)*S3RTF(AKIS)))) 
BETA1    *    RBETAW130./3. 1415927 

12  CONTINUE 

OXSOR    *    (RSTARM l.+EFF))*( (PHI »PHI )• (1 .-RSTAR) »(PHIR- ( VELRAT* 
VELRAT*C0SF(AL=A1)*C0SF(ALFA1)/(PSI*PSI)))-(2.*PHIR»PHI 

2  *S  I\JF  (  ALFA  1 )  •SQ*  TF  (1. -RSTAR) /SQRTF(AKIS))  +  (  (PHI  R*R  ADR  AT* 

3  RADRAT-1. )/AKIS )  ) 
EXX    =    SQRTF(XSQR) 

RBETA2=ATANF(-PSI  «EXX/ (PHI*VELRAT*COSF( ALFAl ) »SQRTF( l.-RSTAR)) ) 

BETA2    *    RBETA2*130./3. 1415927 

OBETA    *    DETA1-BETA2 

PSIA    *    0.99    -    .030228*DBETA    -    U. 97/ ( 1 80. -OBETA) 

DPSI    *    ABSF(PSI    -PSIA) 

IF    (.0000001-OPSI )     10,11,11 

10  PSI    =    PSIA 
K    *    K  +  l 

IF     (K-1000)     12,12 ,45 

11  PSI    »    PSIA 

OETA    »    (2./SQRTFUKIS) )*( ( PHI*SI NF ( ALFAl ) *SQRTF( 1 .-RSTAR ) )- 

1  (RADRAT»RADRAT/SORTF(AKIS) )♦ ( PSI»EXX*RADRAT ) ) 
IF    (ETA-ETAA)     13, 14,14 

13  DETA    *    ABSF(ETA-ETAA) 

IF    (.0000001-DETA)     15,16,16 

14  CONTINUE 
ETAA  -    ETA 
DKIS  *  (-.1)**J 
AKIS  »  AKIS+DKIS 
GO  TO  17 

15  J  »  J*l 

IF  (J-5)  14,14,43 
43  PRINT  44 

GO  TO  50 
•45  PRINT  46 

GO  TO  50 

16  CONTINUE 

OAKE  »AKIS*( (PHI*>HI»( 1 .-RSTAR )»COSF ( ALFAl ) #COSF( AL FA1 ) «VELRAT* 
VELRAT)-M(PSI»PSI)»XSQR)-(2.*PSI»EXX#RADRAT/SQRTF(AKIS) ) 

2  ♦  ( (RAORAT«RADRAT)/AKIS) ) 

VELM1  *  PHI «SQRTr (1 .-RSTAR) *SQRTF( AKIS )»COSF( ALFAl) 
0ALFA2  *  ATANFU  RADRAT/SQRTFt  AKIS)  )-(PSI»EXX)  )/( PH I» VELRAT* 

1  COSFtALFAl )»SQRTF( l.-RSTAR))) 
ALFA2  *  ALFA2*183./3. 1415927 

0WS3R2  =  (PSI*PSI»AKIS)*(  (  (1.+EFF)*RSTAR)-MPHIR*PHI*PHI» 
( l.-RSTAR) )-(2.*PHIR*PHI*SQRTF( 1 .-RSTAR ) *S INF (ALFAl ))/ 

2  (SORTF(AKIS)  )  *■   (  PHI  R+R  A0RAT«RA0RAf-l  .) /AKIS  ) 

QWSQR1  -  1  .MAKlS»PHl*PHt»(  l.-RSTAR)  )*-2.*SQRTFlAKiS  )»PHJ« 
1    SQRTF( l.-RSTA* )*SINF( ALFAl ) 

RELRAT  *  SQRTF(W5QR2)/SQRTF(WSQR1) 
OPRINT  4 1,NALFA1,M PHI E, RADRAT, VELRAT, PHI , PHIR, PSI ,^ STAR, AKIS, 
1ETA,AKE,VELM1,REIRAT,BETA1,BETA2.ALFA2 

L  «  1 

■  • 
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19  CONTINUE 

ETAR  »  ETA-.005 
2U  AKIS«AKIS*.l 
K=l 
ORBETA1  =  ATANF(TANF(ALFA1)  -  ((  1.0)  /  ( PHI»COSF ( ALFA  1 ) •SQRTF 
1   (1.  -  RSTAR)*S3RTF( AKIS)) ) ) 
BETA1  =  RBETAl»130./3. 1415927 

22  CONTINUE 

OXSQR  =  (RSTARft(1.+EFF))*(  ( PHI*PHI )•( 1 .-RSTAR )•( PHI  R-( VELRAT* 

1  VELRAT«C0SF(ALFA1)«C0SFCALFA1 )/(PSI»PSI)) )-(2.*PHIR»PHI    . 

2  «S INF (ALFA1)»SQ*TF(1.-RSTAR) /SQRTF ( AKIS ))+((PHIR*RADRAT« 

3  RADRAT-1. )/AKIS ) ) 
EXX    =    SQRTF(XSQR) 

RBETA2*ATANF(-PSI«EXX/(PHI»VELRAT«C0SF(ALFA1 ) ♦SQRTF (1 .-RSTAR ) ) ) 
BETA2    *    RBETA2M30./3. 1415927 

DBETA  =  BETA1-BETA2 

PSIA    *    0.99    -    .0D0228«DBETA   -    4.97/ < 180. -DBETA) 

DPSI  =  ABSF(PSI  -PSIA) 

IF  (.0000001-DPSI )  21,25,25 
21  PSI=?SIA 

KsK  +  1 

IF  (K-1000)  22,22,U5 
25  PSI=PSIA 

OETA    =    (2./SQRTFUKIS) )•( ( PHI*SI NF ( ALFAl ) *SQRTF ( 1 .-RSTAR) )- 
1        (RADRAT«RADRAT/SQRTF(AKIS)  ) -»■(  PS  I«EXX«RADRAT  ) ) 

IF     (ETA-ETAR)    23,24,24 

23  CONTINUE 
M=1 

55    AKISsAKIS-.OO'l 
M*M*1 
IF    (1000-M)    51,35,35 

35  CONTINUE 

0RBETA1  *  ATANF(T*NF(ALFA1)  -  ((  1.0)  /  ( PHI*COSF ( ALFAl ) »SQRTF 
1   (1.  -  RSTAR)*S3RTF( AKIS)) )) 
BETA1  =  RBETAl*130./3. 1415927 
OXSQR  =  (RSTAR«( 1. +  EFF) )♦( (PHI«PHI )*(1 .-RSTAR) ♦( PHI R-( VELRAT* 

1  VELRAT»COSF( ALFAl ) «C0SF ( ALFAl ) / ( PS  I *PS I ) ) )- (2. •»HIR«PHI 

2  •SINF(ALFA1)*SU*TF(  1  .-RSTAR)  /SQRTF  (  AKI  S)  )♦(<  PHI  R*RADRAT* 

3  RAORAT-1 . )/AKIS ) ) 
EXX  s  SQRTF(XSQR) 

RBETA2=ATANF(-PSI«EXX/(PHI*VELRAT«C0SF(ALFA1)«SQRTF( 1. -RSTAR)) ) 
BETA2  -  RBETA2»130./3. 1415927 

DBETA  =  BETA1-BETA2 

PSIA  =  0.99  -  .0D0228»DBETA  -  4.97/ ( 1 80. -DBETA) 

DPSI  =  ABSF(PSI  -PSIA) 

IF  (.0000001-DPSI)  36,38,38 

36  PSI=PSIA 
GO  TO '35 

38  PSI=PSIA 

OETA  =  (2./SQRTFUKIS)  )*(  ( PHI »SI NF (ALFAl  ) *SQRT F( 1 .-RSTAR )  )- 
1   ( RADR AT »RADRAT/ SQRTF ( AKIS) ) ♦ ( PSI *EXX»RADRAT ) ) 
IF  (ETA-ETAR)  55,26,26 

260 AKE  -AKIS«( (PHI»*HI»( 1  .-RSTAR ) «C0SF( ALFAl ) ♦C0SF( AL FA1 ) ♦VELRAT* 

1  VELRAT )♦(( PSI*PSI) »XSQR)-(2.»PSI*EXX«RADRAT/SQRTF( AKIS) ) 

2  +T(RADRAT*RADRAT)/AKIS)) 

VELM1  =  PHI*SQRTF ( 1 .-RSTAR) *SQRTF ( AKI S )*C0SF( ALFAl  ) 
0ALFA2  *  ATANF(  (  U ADRAT/SQRTF (  AKI  S)  )-(  PSI  »EXX)  )  /  (  P-l  I*VELRAT* 

1  C0SFULFA1  )»SQRTF(1. -RSTAR) ) )  . 
ALFA2  =  ALFA2*183./3. 1415927 

0WSQR2  -  (PSI*PSI»AKIS)*(  ( ( 1  . +EFF ) «RSTAR) ♦ ( PHI R*PHI *PHI* 

(1. -RSTAR) )-(2.«PHIR»PHI«SQRTF(1.-RSTAR)*SINF(ALFAl ))/  , 

2  (SQRTF (AKIS))* ( PHIR*RADRAT»RADRAT-1 . ) /AKIS ) 

0WSQR1  =  1.+(AKIS»PHI*PHI*(1.-RSTAR))-2.*SQRTF(AKIS)*PHI* 
1    SQRTF( l.-RSTA* )»SINF(ALFA1 ) 

RELRAT  =  SQRTF(WSQR2)/SQRTF(WSQR1) 
0  PR  I  NT  41,NALFA1 ,^ PH IE, RADRAT, VELRAT , PHI ,PHIR, PSI, *STAR, AKIS, 
1ETA,AKE,VELM1 ,REL RAT , BETA1 , BETA2, ALFA2 

L=L  +  1 

IFU-3)  19,42,42 
42  PRINT  47 

-.  £Q  T0  §2 
51  PRINT  52 

GO  TO  50 
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U00F0RMAT(/n8Hl  Al    PE    R2/R1    VM2/VM1    PHI    PHIR     PS  I     R« 

KIS       ETA      KE    VM1/U1   W2/W1    BETA1    BETA2   ALFA2   /) 
410F0RMAT  (/I3,I6.F7.1,F9.1tF8.2,F7.2,2F7.3,F8.3tF9.5t2F7.3. 

1  Ffi  ^  ?F3  ?  FT  ?  ) 
44  FORMAT  (/35H  FIVE  ITERATIONS  OF  AKIS  COMPLETED) 

46  FORMAT  (/3UH  100D  ITERATIONS  OF  PSI  COMPLETED  ) 
52  FORMAT  (35H  1000  ITERATI ONS  .OF  AKIS  COMPLETED    ) 

47  FORMAT  (///) 
4G  FORMAT  (1H1) 

50  CONTINUE  '  • 

END 
END 
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PROGRAM  FANN03 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


A 

P/P 

DAB 

DAN 

EE 

EN 

n 

EXP1 

EXP2 

EXP3 

EXP4 

GAM 

v 

PHI 

<t> 

PH12 

<fi2 

RA 

RB 

R 

ZETA 

1 

ZETAP 

f, 

Pressure  ratio 

Name  for  the  term  (p/P) 

expl 
Name  for  the  term  (p/P) 

Reheat  factor  defined  by  (  f    -  }  )/(l  -  7  ) 

Polytropic  exponent.   See  Appendix  I 

Exponent  #1  »  (2/n) 

Exponent  #2  •  (n+l)/n 

Exponent  #3  -  (1  -  f  p)  (  tf  -  1)/^ 

Exponent  #4  -  (tf  -  1)/^ 

Specific  heat  ratio 

Flow  function 

Flow  function  squared.   Intermediate  step  in  finding 

6x0  3 
Name  for  the  term  (p/P) 

6XdA 

Name  for  the  term  (p/P) 

Computed  loss  coefficient  to  be  used  for 
determining   p.   R  -  (RA-RB)/(1-RB) 

Loss  coefficient  ■  1  -  y 

Polytropic  loss  coefficient  defined  in  Appendix  I 
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PROGRAM    FANN03 

DIMENSION    PHM  11)  ,ZETA<11  ),ZETAP(  1 1  ),EE(  11) 
DO    22    1    *    125, JUD 
XI    *    I 

GAM    »    XI/IOO. 
100F0RMAT    (26H1FLOW    FUNCTION    FOR    GAMMA    *    F5.2// 
19H    PRESSURE    ) 
3100FORMAT    I30H1    ZETA    POLYTROPIC    FOR    GAMMA    «    F5.2// 

1 9H    PRESSURE    ) 
2100F3RMAT    (32H1    REHEAT    FACTOR    »E«    FOR    GAMMA    »    F5.2// 
19H    PRESSURE    ) 
WRITE    OUTPUT    TAPE    3,10.GAM      . 
WRITE    OUTPUT    TAPE    M, 310, GAM 
WRITE    OUTPUT    TAPE    5, 210, GAM 
DO    11    M    *    1,11 

11  ZETA(M)    «    -.025*.025«FLOATF(M) 

12  FORMAT    (18H       RATIO  ZETA    *    11F9.3/) 
WRITE    OUTPUT    TAPE    3, 12, ZETA 

WRITE    OUTPUT    TAPE    U, 12, ZETA 
WRITE    OUTPUT    TAPE    5, 12, ZETA    . 
A    «    1,02 

13  CONTINUE 

00    111    N    *    It  11 

ZETA(N)    *    -.025*.025»FL0ATF(N) 

ZETAP(N)    »    -.001 
102    ZETAP(N)    «    ♦.OOHZETAP(N) 

EXP3    *    (1.-ZETAP(  N))MGAM-1.)/GAM 

EXPU    *    (GAM-1.)/3AM 

RA    *    1./(A»*EXP3) 

RB    «    1./(A#«EXPU) 

R    *    (RA-RB)/( l.-RB) 

DIFF    »    ZfcTA(N)-R 

IF    (OIFF)    101,103,102 
101    ZETAP(N)    *    -.OOODUZETAP(N) 

EXP3    *    (l.-ZETAP(N))«(GAM-l.)/GAM 

EXPU    -    (GAM-D/3AM 

RA    *    l./(A»»EXP3) 

RB    «    l./(A«»EXPU) 

R    *    (RA-RB)/( l.-*B) 

OIFF    »    ZlTA(N)-R 

IF  (DIFF)  101,103,100 
100  E^  »  BAM/C l.*ZET*P(N)»(GAM-1.)) 

EXP1  «  2. /EN 

EXP2  *  (EN>1. )/EM 

DAN  «  1,/(A**EXP1  ) 

DAB  *  l./(A**EXP2) 

PHI2  »  2.*(GAM/OAM-l.))«(DAN-DAB) 

PHKN)  »  SQRTFCP-II2) 
111  EC(N)  *  (ZETAP(N)-ZETA(N))/M.-ZETA(N)) 

WRITE  OUTPUT  TAPE  3, 15, A, PHI 

WRITE  OUTPUT  TAPE  U,15,A,ZETAP 

WRITE  OUTPUT  TAPE  5,15.A,EE 
15  FORMAT  (F6.2.13X.11F9.5) 

IF  (A-1.99)  17,19,19 
17  A  »  A*. 02 

GO  TO  13 

19  IF  (5.9-A)  22,20,20 

20  A  *  A*. 20 

21  GO  TO  13 

22  CONTINUE   . 
400  FORMAT  (8FT2.5) 

END 
ENO 
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PROGRAM  FANN04 


FORTRAN  NAMES,  EQUIVALENT  SYMBOLS  AND  DEFINITIONS 


A 

CPHI 

DAB 

DAN 

DIFF 

EN 

n 

EXP1 

EXP4 

EXP5 

EXP6 

GAM 

* 

PH 

(P/P) 

PHI 

<P 

PHI  2 

<? 

PHIC 

<tc 

PHICR 

rt. 

PPCR 

(p/p) 

PR 

p/p 

R 

RA 

RB 

cr 


cr 


Interim  critical  pressure  ratio  while  computing 
Compute  <P   ,  subroutine  name 

Name  for  the  term  (p/P) 

exp5 
Name  for  the  term  (p/P) 

Difference  between  given  and  computed  loss 
coefficient 

Polytropic  exponent 

Exponent  #1  -  {t   -  1)/-? 

Exponent  #4  -  (1  -  ^p)  (  tf  -  l)/T 

Exponent  #5  ■  2/n 

Exponent  #6  ■  (n+l)/n 

Specific  heat  ratio 

Critical  pressure  ratio  name  for  storing  away 
and  printing 

Flow  function  assumed  for  comparison 

Flow  function,  squared 

Flow  function  computed  in  subroutine  based  on 
tf  and  *5 

Choked  flow  value  of  flow  function 

Choked  flow  value  of  pressure  ratio 

Pressure  ratio 

Computed  loss  coefficient 

Name  for  the  term  (p/P)exp 

Name  for  the  term  (p/P) 
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' 


ZETA 
ZETAP 


Polytropic  loss  coefficient  name  for  storing  away 
and  printing 

Loss  coefficient 

Polytropic  loss  coefficient 
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PROGRAM    FANN04 

DIMENSION       ZETAdl  )  , PR( 1 1 ) ,PH(1 1 ) , Z ( 1 1 ) 

PRINT  UOO 

DO  U99  J=l,ll 

499  ZETA( J)*-.025+.025»FLOATF(J)  . 
PRINT  U01,ZETA 

DO  U90  I  *  125,1*0 
XI  =  I 

GAM*  XI/100. 
DO  498  K-1,11 
A  =  1.8 
PHI=0.2 
U97  CONTINUE 

CALL  CPHI  (GAM, Z5TA(K), A. PHIC, PHICR, PPCR, ZETAP) 
IF  (PHIC-PHI)  49i»,493,495 

495  CONTINUE 
PHI  -    PHIC 
A  *  PPCR 
GO  TO  497 

49U  DIFF  *  ABSF(PHIC-PHI) 

IF  (OIFF-. 0000001 )  493,493,496 

496  CONTINUE 
PHI  =  PHIC 
A  *  PPCR 
GO  TO  497 

4,93  CONTINUE 

PR(K)  =  PPCR 

PH(K)  *  PHICR 

Z(K)  *  ZETAP 
4  98  CONTINUE 

IF(  1-133)  492,491  ,492 

491  PRINT  400 
PRINT  401 ,ZETA 

492  PRINT  402, GAM 
PRINT  412, PH 
PRINT  403, PR 
PRINT  404.Z 

4000F3RMAT  (90H1  CHO<ED  FLOW  VALUES  OF  THE  FLOW  FUNCTION,  PRESSURE  RA1 
110  AND  POLYTROPi:  LOSS  COEFFICIENT    //) 

401  FORMAT  (11X7HZET4  *  11F9.3//) 

402  FORMAT  ( 7H  GAMMA*  F4.2) 

412  FORMAT  ( 1 8H  PHI  *   11F9.5/) 

403  FORMAT  ( 1 8H         (P0/P1C  *  11F9.5/) 

404  FORMAT  ( 1 8H  ZETAP  »  11F9.5//) 
490  CONTINUE 

END 

SUBROUTINE  CPHI  ( GAM, ZETA, A, PHIC, PHICR, PPCR, ZETAP) 

ZETAP  =  -.001 

500  ZETAP=0.001+ZETA> 

EXP4M  1. -  ZETAP  )*(G  AM- 1.)  /GAM 
EXP1    =    (GAM-1 .J/SAM 
RA    =    l./(A»»EXP4) 
RB    -     l./(A»«EXPl) 
R    =    ( RA-RB)/{ l.-*B) 
IF(ZETA-R)    501,532,500 

501  ZETAP*-. 00001+ZETAP 

EXP4=( l.-ZETAPU(  GAM-1.) /GAM 
EXP1  =  (GAM-1. J/JAM 
RA  "     l./(A**EXP4) 
R3  *  l./(A»«EXPl) 
R  *  (RA-RB)/( l.-*B) 
IF(ZETA-R)  501,532,502 

502  EN=GAM/( 1 .*ZETAP» (GAM-1.  )) 
EXP5  *  2. /EN 

EXP6  =  (EN+1. )/EM 

DAN  *  l./(A»»EXPS ) 

DAB  -    l./(A»»EXP6) 

PHI2  =  2.«(GAM/(3AM-1.))»(DAN-DAB) 

PHIC  -    SQRTF(PHI2) 
0PHICR*((2./(EN*1.  ))»«(1./(EN-1.)))«S0RTF((2./EXP1)«((EN-1.)/ 
I  (EN+LH) 

PPCR  *  ((EN*1.)/2.)«»(EN/(EN-1.)) 

ENO 
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PROGRAM  AFTER 

AME  =1.37 

NALFA  -    75 

ALFA    -    FLOATF(NALFA)»3.1415927/180« 

GAMMA    «    1.37 

EXP  *  (GAMMA  -1.) /GAMMA 

EE  ■  (2. /((GAMMA-1. )»AME»AME)) 

PRINT  409 

PRINT    410,AME,NALFA,GAMMA 

PRINT  404- 

DO  50C  I  '    1,26 

PRAT  *  1.02  -  .02»FL0ATF(L) 

DALFA  =  0.0 

M  =  0 

100  M  -    M*1 

101  CALFA  *  (-.1  )»»(M-1)+0ALFA 
DALFAR  =  DALFA*3. 1415927/180. 
COSRAT  =  COSF(ALFA)./CCSF(ALFA-OALFAR) 

AA  *  GAMMA»AME*AME*C0SF(DALFAR)-COSRAT«(GAMMA-l.)»AME»AME/2. 
BB  *  -GAMMA»AME»AME-1. 

CC    =    COSRAT«( l.*(  (GAMMA-1 .  )»AME*AME/2 . ) ) 
DD    =    EB»5B-4.»AA«CC 
IF    (DC)    101,200,200 
200    VELRAT1    =    (-BB+SQRTF ( DD ) ) / ( 2. »AA) 

0PRAT1    =    ( COSRAT* (  l.-M  (GAMMA-1. )»AME»AME/2. )•< l.-VE LRAT 1*VELRAT 1 )  )  ) 
1       /VELRAT1 
2A    *    1.-(VELRAT1*VELRAT1)/(1.*EE-EE*(PRAT1««EXP)) 
DIFF    =    ABSF(PRATI-PRAT)     u 
IF    (DIFF-. 0001)300,300.102 

102  IF    (PRAT1-PRAT)    103,300,101 

103  M    *    M*1 

104  DALFA    =     (-. 1)*»(M-1 )+DALFA 
DALFAR    =    DALFA»3. 1U15927/180. 
COSRAT    *    COSF(ALFA)/CCSF(ALFA-DALFAR) 

AA    =    GAMMA*AME»AME»C0SF(CALFAR)-C0SRAT*(GAMMA-1.)»AME*AME/2. 

BB    =    -GAMMA«AME*AME-1 . 

CC    -    COSRATM  1.-M  (GAMMA-1. )»AME*AME/2.)) 

DD    =    eB«BB-U.*AA»CC 

IF     (DC)     100,105,105 

105  VELRA.T1    =     (-PB+SQRTF  (  DD)  )  /  (2.  *AA) 

0PRAT1    *    (COSRAT«(  1 .♦(( GAMMA-1 .) •AME«AME/2. )•( 1 .-VELRAT1 «VELR AT  1 )) ) 
1       /VELRAT1 

ZA    *    1.-( VELRAT1»VELRAT1)/(1.*EE-EE«(PRAT1««EXP) ) 

DIFF    *    ABSF(PRATl-PRAT) 

IF    (DIFF-. 0001)    300,300,106 

106  IF    (PRAT-PRAT1)    100,300,104 

300    PRINT    U12,PRAT,PRAT1,VELRAT1, DALFA, ZA,AMO,AMA 

GO    TO    500 
401     FORMAT    (20H    B2-4AC    IS    NEGATIVE         /) 
U020F0RMAT    (60H  ABC  CCSRAT  GAMMA 

1    AME  ) 

U03    FORMAT    (6F10.6/) 

404  FORMAT    (U6H         PRAT  PRAT1  VELRAT1  OALFA  ZA       /) 

405  FORMAT    (/E20.10) 

406  FORMAT  ( F7.4, 4F 10 .4 ) 

409  FORMAT  (25H1     M ACH   ALFA      GAMMA     /) 

410  FORMAT  ( F  10.3, 1  5, F1 0. 3  /) 

411  FORMAT  (66H       PRAT      PRAT1    VELRAT1      DALFA         ZA 
1    MD     MA    /) 

412  FORMAT  (7F10.4) 
500  CONTINUE 

END 
END 
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Appendix  IV-   Critical  Mach  Number  of  Flow  Undergoing  Polytropic 
Expansion. 


The  critical  pressure  ratio  for  a  flow  undergoing  a  polytropic 
expansion,  as  developed  in  Appendix  I,  is 


P_ 


z 


1*, 


(1,13) 


To  find  the  maximum  velocity  corresponding  to  this  pressure 
ratio  it  must  be  recalled  that,  for  a  polytropic  process 


-M. 


-4>  4f~      =  constant 


H-fT- 1 


where  C  is  a  constant 


t-m        sn 


constant 


(IV, 1) 


(IV, 2) 


Therefore 


r_ 


P, 


Kff  (IV, 3) 

Using  the  T-s  diagram  of  Fig.  24  as  a  reference,  it  can  be  seen  that 


Z-t  =  r0 
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and 
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(IV, 4) 


For  the  critical  velocity  then,  using  (1,13) 
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which  reduces    to 
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The  acoustic  velocity  is 

A*. 
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which  reduces   to 
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Therefore 
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Also  developed  in  Appendix  I  was  the  relation  for  the  polytropic 
exponent  n,  in  terms  of  the  polytropic  loss  coefficient,  as 


and 


M  /-  v 

'  is*     ss 


'&L 


I  *■  ff(r  -0 


For  example,  for  3*-  1.37  and  *f     -  0.10,  M   ■  .931, 

and  for   ^   -  0.20,  M   -  .854. 
5  P         « 
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(1,7) 


Making  this  substitution,  and  collecting  terms 

'f  (IV, 8) 


(IV, 9) 


